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20.  Abstract  (Continued) 

survival  based  on  a simulated  missile-target  engagement,  including  the 
dynamics  anti  principal  characteristics  of  the  aircraft  and  the  missile. 
The  program  can  also  simulate  the  deployment  of  decoys  such  as  flares 
or  pyrophorlcs. 

The  utility  of  the  program  lies  in  that  it  can  provide  guidelines 
during  aircraft  configuration  studies,  assess  effects  of  design  changes  on 
aircraft  survivability,  and  permits  tradeoff  studies  to  be  made  between 
various  CMs  such  as  suppression,  shielding  and  flare  deployment. 

The  programs  are  operational  on  the  CDC  6600  digital  computer 
at  Wrlght-Patterson  Air  Force  Base,  Ohio. 
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INTRODUCTION 


Thl*  document  constitutes  the  final  report  on  Analytical  Methodology 
for  Evaluation  of  Payoffs  for  Infrared  Countermeasures  and  Suppression 
(EPICS).  Contract  No.  F33615-75-C-4076.  The  prime  objective  of  this 
study  was  to  develop  a methodology  or  analytical  tool  for  rapidly  and  effi- 
ciently assessing  the  Impact  of  Infrared  suppression  techniques  on  aircraft 
survivability.  Specifically,  the  intended  purpose  of  the  methodology  is  to 
provide  a capability  to  analytically  predict  the  effectiveness  of  aircraft  design 
changes  (primarily  those  related  to  infrared  signature)  on  the  probability  of 
sirersft  survival  in  a specified  infrared  threat  environment. 

The  development  of  an  analytical  tool  that  meets  the  above  objectives 
was  achieved.  This  tool  consists  of  two  complementary  digital  computer 
programs:  (1)  an  Infrared  target  signature  model  (ASDIR  II)  and  (2)  a missile/ 
target/countermeasures  (M/T/CM)*  model.  A third  program  SPRINT  (a 
subroutine)  provides  the  interface  between  the  two  models.  All  three  pro- 
grams are  fully  operational  on  the  CDC  6600  computer  system.  This  total 

methodology  system  is  designated  as  EPICS. 

The  first  program,  ASDIR  II,  was  primarily  developed  by  the  Air  Force. 
Hughes  modified  it  and  made  it  operational.**  It  is  documented  under  a 
separate  cover.***  The  Inputs  to  ASDIR  11  are  engine  specification  data  (gas 
dynamics  or  measured  plume  data  to  determine  the  engine  exhaust  plume 
radiation)  and  engine  hot  metal  parts  in  terms  of  temperature  and  radiating 

♦The  baseline  for  the  M/T/CM  program  was  developed  by  Hughe* > under 
"IRCM  Simulation  Study",  Contract  No.  F33615-C-74-1680  for  AFAL. 

**This  task  constituted  a deviation  from  the  Statement  of  Work,  Originally, 
a Hughes' -developed  target  signature  program.  IRSIG,  was  to  be  used. 
However,  on  Program  Monitors  instructions,  ASDIR  II  was  used  instead. 

♦♦♦Stone,  Charles  W. , Capt. , USAF  and  Tate,  Stanley,  ASDIR  II  (Vol.  I,  II, 
and  III).  Deputy  for  Development  Planning,  Aeronautical  Systems  Division, 
Wrlght-Patterson  Air  Force  Base,  Ohio,  No.  ASD/XR-TR-75- 1,  January 
1975. 
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area  ah  a function  of  aspect.  Similarly,  the  contribution  to  the  total  IR 
signature  due  to  skin  aerodynamic  heating  is  an  input  in  terms  of  temperature 
and  radiating  area  for  as  many  as  twenty  skin  surfaces.  Finally,  the  viewing 
geometry  — target  and  observer  altitudes,  aspect  angle,  and  slant  range  is 
an  input  to  the  AS  DIR  11  program. 

The  outputs  of  this  code  are  in  the  form  of  polar  data  for  the  source 
spectral  radiant  intensity,  J^,  integrated  over  the  missile  response  bar.d 
and  the  apparent  spectral  radiant  intensity,  J^r^,  *ls°  integrated  over  the 
missile  band,  but  at  the  point  of  a remute  observer.  These  data  then  serve 
as  input  to  the  M/T/CM  program. 

The  second  major  element  of  EPICS  is  the  M/T/CM.  This  program  is 
described  in  detail  in  this  report.  The  program  is  a generic  five-degrees- 
of-freedom  dynamic  simulation  of  the  total  missile/target  encounter  in  a 
countermeasures  environment.  The  prime  countermeasures  technique  that 
can  be  evaluated  using  this  program  are  IR  signature  reduction  through 
suppression  or  shielding,  and  active  decoys  such  as  flares  or  pyrophorics. 
The  output  of  the  program  is  a probability  of  target  survival  (Pg)  under  a 
varied  set  of  launch  conditions  and  for  various  IR  missile  threats.  The  Pg 
is  defined  by 


number  of  misses 

total  number  of  launch  cases 


As  indicated  in  a preceding  footnote,  the  baseline  subroutines  for  the 
M/T  /CM  program  were  developed  by  Hughes  under  an  earlier  Air  Force 
study  program,  however,  in  the  present  study  contract  this  baseline  pro- 
gram was  considerably  expanded  and  improved.  In  addition,  the  program 
was  modified  to  accept  inputs  from  the  ASDIR  II  program  with  the  aid  of  a 
subroutine  called  SPKINT.  The  total  program  was  compiled  on  the  CDC  6600 
computer  system  (it  was  originally  written  for  the  SIGMA  5 Computer).  The 
major  changes  to  M/T/CM  include: 

Modularization  of  the  program 

Addition  of  flare  control  options  (function  of  range  and  time-to-go) 


x 
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AUdition  of  superelevation  angle  subroutine  for  ground-to-air  missiles 

Addition  of  launch  mode  selection  option  in  azimuth  and  elevation 

Reduction  of  program  execution  time 

Sample  runs  using  all  three  programs  were  also  made;  see  Section  10. 

In  this  program,  a specific  target  is  represented  by  Us  physical  char- 
acterietice,  its  dynamic  parameters,  and  its  Infrared  signature.  The  physi- 
cal characteristics  include  the  size  of  the  aircraft,  its  wing  span  and  the 
engine  location.  The  dynamic  parameters  delude  initial  velocity,  accelera- 
tions, and  maneuvers.  The  infrared  target  signature  is  represented  by  the 
effective  apparent  radiant  intensity.  Jt^.  ns  a function  of  range  and  aspect 
angle,  and  is  calculated  by  ASD1R  U.  Interpolation  on  range  and  aspect 
angle  between  (Jt)*x  data  points  provide  the  appropriate  values  during  the 
simulated  flight  which  in  turn  ave  used  to  calculate  the  effective  irradiance 
at  the  missile  seeker. 

Threat  missiles  are  represented  by  a number  of  parameters  divide  in  o 
six  categories;  *eeker,  signal  processing,  guidance,  aerodynamics,  propul- 
sion, and  physical  characteristics.  Currently,  12  missiles,  25  aircraft,  and 
4 flare  types  have  been  defined  and  are  part  of  the  simulation  file. 

The  M/T/CM  program  has  been  validated  by  comparing  simulated  engage- 
ment. o*  captive  missiles  being  decoyed  by  flares,  with  flight  test  results 
(using  same  flight  conditions  and  the  same  missiles)  conducted  by  the  Naval 

Weapons  Center,  China  Lake,  California. 

As  mentioned  above,  the  interface  between  ihe  AS  DIR  II  signature  pro- 
gram and  the  M/T/CM  program  is  provided  by  an  auxiliary  spectral  integra- 
tion subroutine,  called  SPKINT.  This  routine  integrates  the  apparent  spectral 
radiant  intensity  values.  Vx . over  any  specified  spectral  interval  to  obtain 
effective  radiant  intensity  (Jt)aX.  In  general,  the  integration  is  performe 
for  spectral  Intervals  corresponding  to  the  spectral  bandpasses  of  the  12  mis- 
siles on  file.  The  SPKINT  subroutine  is  described  in  Section  9. 

In  summary,  the  EPICS  methodology  provides  a tool  to  assess  the  impact 
that  aircraft  dosign  has  on  the  aircraft  survivability  in  an  infrared  missile 
threat  environment,  evaluating  design  changes  and  conducting  tradeoff  studies 
during  preliminary  design  and  determining  aircraft  survivability  in  a flare 
countermeasures  environment. 

xi 
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Figure  1 ahowa  a flow  diagram  for  the  aimulation  program  executive 
operation  that  calla  all  subroutines  The  program  ia  broken  into  eight  major 
areaa  with  each  area  being  subsequently  discussed  in  Sections  1 through  8 of 
this  report. 

Section  1 deala  with  the  creation  of  files  on  which  the  necessary  constants 
for  the  missile,  target,  and  flare  are  stored.  Long  lists  of  input  data  can  be 
eliminated,  by  creating  files  for  each  missile,  target,  and  flare  to  be 
evaluated  and  the  required  constants  can  be  specified  by  simply  referring  to 
a file  name. 

In  Section  2,  the  launch  geometry  variables,  the  flare  deployment 
strategy,  tire  aircraft  maneuver  option,  and  all  other  program  options  are 
set.  In  this  section,  all  program  variables  are  initialized. 

Section  3 of  the  program  updates  the  position,  velocity,  and  acceleration 
of  the  missile,  target,  and  flare(s)  with  respect  to  inertial  coordinates. 

In  Section  4,  relative  ranges,  range  rates,  angles,  and  angular  rates 
are  determined  between  the  missile  and  the  target  as  well  lb  the  missile  and 
flare(s). 

The  irradiance  at  the  missile  dome  from  the  target  and  the  flare(s)  is 
computed  in  Section  5. 

In  Section  6,  the  aimpoint  location  is  determined  based  on  the  irradiance 
levels  of  the  sources  in  its  field  of  view  (FOV)  and  the  type  of  signal  process- 
ing in  the  missile.  This  aimpoint  location  is  then  fed  back  into  the  missile 
dynamics  (Section  3)  through  missile  guidance. 

When  the  simulation  prog  ran.  goes  into  an  abort  mode,  the  point  of 
closest  approach  of  the  missile  to  the  target  is  determined.  The  details  are 
given  in  Section  7. 

Section  8 describes  how  the  probability  of  hit  is  determined  based  on  the 
closest  approach  distance,  aircraft  cUmensicnB,  type,  and  lethality  of  the 
missile  warhead. 

The  geometry  for  the  missile  and  target  encounter  is  shown  in  Figure  2. 
This  simulation  is  based  on  a two -plane  geometry,  because  a missile 
essentially  processes  its  target  position  and  rate  information  and  provides 
guidance  commands  in  two  separate  planes  --  horizontal  and  vertical. 

Coupling  between  the  two  planes  is  accomplished  by  the  range  and  velocity 
variables. 
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Figure  1.  Simulation  program  executive  operation  flow  diagram 
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Figure  2.  Missile  and  target  geometry 
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The  coordinate  ayatem  in  which  miaailei  target,  and  flare  poaitiona  are 
calculated  ia  ahown  in  Figure  2.  Thia  ayatem  haa  aa  ita  crigin  a point  on 
the  ground  directly  below  the  launch  point  of  the  miaaile.  The  Z axis  ia 
parallel  to  gravity  and  poaitive  up.  Therefore  the  initial  poaition  of  the 
miaaile  ia  given  by  (O,  O,  Z^),  where  Z ^ ia  the  launch  altitude.  The  X 
axia  ia  perpendicular  to  gravity  and  oriented  auch  that  the  initial  poaition  of 
the  target  ia  in  the  X-Z  plane.  The  Y-axia  completea  the  orthognal, 
right-handed  coordinate  ayatem.  The  initial  poaition  of  the  target  ia  given 
by  (XT»  O,  Z^.)  with  X^.  being  the  horizontal  range  between  the  miaaile  ant1, 
target  at  launch,  and  Zj  being  the  target  altitude.  From  thia  definition,  the 
line-of- eight  (LOS)  ia  in  the  X-Z  plane  from  the  miaaile  to  target  at  launch. 

The  target  aapect  relative  to  the  miaaile  ia  aet  by  the  target  velocity  vector. 

The  atate  vectora  liated  in  Table  1 represent  the  X-,  Y - , and  Z-componenta  - 
of  position,  velocity,  and  acceleration  of  the  miaaile  and  target.  The  state 
variables  are  divided  into  two  seta:  one  set  for  the  horizontal  plane  and  one 
set  for  the  vertical  plane. 

Figure  3 shows  the  geometry  that  exists  between  the  missile  and  an 
arbitrary  flare  (the  Xth  flare).  Again,  as  in  the  case  of  the  miaaile  and  tar- 
get geometry,  the  problem  ia  divided  into  two  planea,  horizontal  and  vertical. 

The  atate  vector  for  the  Xth  flare  ia  listed  in  Table  2.  Aa  in  the  case  of 
the  missile  and  target  atate  vector  system,  the  first  five  components  repre- 
sent the  position,  velocity,  and  acceleration  of  the  miaaile;  while  the  next 
four  components  represent  the  position  and  velocity  of  the  X flare.  It  is 
assumed  that  the  flare  haa  no  thrusting  device,  and  therefore  no  thrust 
acceleration  terms  are  possible.  As  before,  the  atate  variables  are  divided 
into  the  sets  --  one  for  each  plane.  The  executive  routine  listing  for  the 
simulation  program  is  presented  in  Table  3.  All  subroutines  in  the  simula- 
tion program  are  called  by  this  routine. 
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Table  1.  Missile  and  target  state  variable  definitions 


Horizontal  Plane 

Vertical  Plane 

X(l) 

= X - Position  of  Mis-ile 

XP(1) 

= X - Position  of  Missile 

X(2) 

= X - Velocity  of  J 

XP(2) 

= X - Velocity  of  Missile 

X(3) 

= Y - Position  of  Missile 

XP(3) 

= Z - Position  of  Missile 

X(4) 

« Y - Velocity  of  Missile 

XP(4) 

= Z - Velocity  of  Missile 

X(5) 

s Normal  Acceleration 
(XY)  of  Missile 

XP(5) 

= Normal  Acceleration 
(XZ)  of  Missile 

X(6) 

s X - Position  of  Target 

XP(6) 

= X - Position  of  Target 

X(7) 

= X - Velocity  of  Target 

XP(7) 

= X - Velocity  of  Target 

X(8) 

= Y - Position  of  Target 

XP(8) 

=.  Z - Position  of  Target 

X(9) 

= Y - Velocity  of  Target 

XP(9) 

X z - Velocity  of  Target 

X(10) 

= Normal  Acceleration 
(XY)  of  Target 

XP(10) 

= Normal  Acceleration 
(XZ)  of  Target 

Table  2.  Miaaile  and  .Kth  flare  atate  variable  definitiona 


Horizontal  Plane  Vertical  Plane 


XF(1,K)  = X - Poaition  of  Miaaile 

XF(2,K)  = X - Velocity  of  Miaaile 

XF(3,K)  = Y - Poaition  of  Miaaile 

XF(4,  K)  = Y - Velocity  of  Miaaile 

XF(5,  K)  = Normal  Acceleration 
(XY)  of  Miaaile 

XF{6,  X)  = X - Poaition  of  Kth  Flare 
XF(7,  K)  * X - Velocity  of  Kth  Flare 
XF(8,  K)  = Y - Position  of  Kth  Flare 
XF(9,K)  = Y - Velocity  of  Kth  Flare 


XFP(1,K)  = X - Poaition  of  Miaaile 

XFP(2,  K)  = X - Velocity  of  Miaaile 

XFP(3,  K)  = 7 - Poaition  of  Missile 

XFP(4,I<)  = Z - Velocity  of  Missile 

XFP(5,K)  = Normal  Acceleration 
(XZ)  of  Missile 

XFP(6,K)  = X - Position  of  Kth  Flare 
XFP(7,  K)  = X - Velocity  of  Kth  Flare 
XFP(8,  K)  = Z - Poaition  of  Kth  Flare 
XFP(9,  K)  = Z - Velocity  of  Kth  Flare 
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Table  3-  Simulation  program  executive  routine 
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(Table  3,  continued) 
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1.  FILE  CREATION 


Information  used  to  create  missile,  target  and  flare  file  data  hat  been 
collected  at  Hughes  over  a long  period.  Pertinent  sources  are  listed  In  the 
references  at  the  end  of  this  section. 

MISSILE  FILE 

The  missile  Information  Is  divided  into  seven  categories:  seeker,  signal 
processing,  guidance  unit,  aerodynamics,  motor,  physical  characteristics, 
and  other  characteristics.  Table  1-1  Is  the  missile  file  subroutine.  See 
Table  1-2  for  sample  file  data. 

Seeker 

The  seeker  Is  defined  as: 

1.  Seeker  look  angle  (deg),  variable  name  SA.  This  Is  the  angle 
about  the  longitudinal  axis  of  the  missile  In  any  plane  that  the 
gyro  Is  free  to  move. 

2.  Max  gyro  rate  (deg/sec.),  variable  name  WLIM. 

3.  FOV  (deg),  variable  name  FOV. 

4.  Gyro  time  lag  (sec. ),  variable  name  TGU. 

Signal  Processing 

The  signal  processing  includes: 

1.  The  detector  bandwidth  of  the  missile,  variable  name  IBNDM.  The 
bandwidth  value  Is  currently  deleted  from  the  missile  page  printout. 

2.  Aimpoint  type  is  selected  on  the  basis  of  geometric  centroid, 
irradiance  centroid,  and  maximum  irradiance. 

The  present  available  missiles  are  divided  into  two  categories:  the  con- 
scan  missile  which  Is  a maximum  Irradiance  tracker,  and  the  spin-scan 
missile  which  Is  an  Irradiance  centroid  tracker.  Geometric  centroid  does 
not  apply  to  the  current  systems 
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Table  1-1.  Mlaaile  file  aubroutine 
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Table  1-2.  Missile  constants  file 
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PHYSICAL  CHARACTERISTICS  SThCR  CHARACTERISTICS 

DIAMETER  (FT!  J.7/12  LIFESPAN  iSECl 

■elQHT  ILRl  >0.3  LAUNCH  VELOCITY  (FT/SCCI 

ha»  FLIRPFR  OEFLECTION  {DEJ>  17.0  R*K  FLIPPER  TRAVEL  COEOl 


» 

Guidance  Unit 

The  guidance  unit  la  composed  of: 

1.  The  navigation  constant,  variable  name  GKO 

2.  No  guidance  period  (sec),  /ariable  wme  TB 

3.  G-limit  (g's),  variable  name  AS 

4.  Missile  time  constant  (sec),  variable  name.  TS. 

Aerodynamics 

A sampling  of  eight  values  has  been  adopted  for  the  parameters  below 
to  cover  the  broad  range  values  in  order  of  convenience  In  handling  data 
Information  for  table  lookups.  (Exceptions  will  La  noted.)  Tne  following 
variables  are  functions  of  Mach  number  table,  variable  name  VMC: 

1.  Chord  force,  variable  name  CDO. 

2.  Base  drag  coefficient,  variable  name  CDB. 

3.  Maximum  normal  force  coefficient,  variable  name  CNT. 

4.  Angle  of  attack  (deg)  (80  samples),  variable  nane  ALPH.  This 
parameter  Is  a function  of  both  mach  number  and  normal  force, 
variable  name  CNA. 

Motor 

The  motor  parameters  are  depicted  as  tables  of  10  samples  related  to 
time.  These  are; 

1.  Time  (sec. ),  variable  name  ST 
?.  Thrust  (lb. ),  variable  name  TL 

3.  Specific  Impulse  (l/sec. ),  variable  name  SIM 

4.  Mou?«  weight  drop  (lb. ),  variable  name  WID. 

See  Table  1-3, 
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Table  1-3.  Motor  characteristics 


Example: 


Time 

0. 

1. 

1.395 

1.405 

a 

n 

5. 

6. 

7. 

8. 

Thrust 

4660. 

4660. 

4660. 

0. 

1 

0. 

0. 

0. 

0. 

Specific  Impulse 

210. 

210. 

210. 

210. 

210. 

210. 

210. 

210. 

Motor  Weight  Drop 

0. 

0. 

0. 

0. 

■ 

ES 

m 

0. 

0. 

0. 

Thrust 

Level 


4660. 


6 2 4 5 8 

Time 


Weight 


134.  2 


X 


.101 


1 4 5 8" 

Time 


Equations 


1.  Thrust  = TL 

2.  Weight  (new)  = Weight  (old)  - Thrust  (DELT/SIM)  - WD(t) 


DELT  = Integration  Step  Size 


Physical  Characteristics 

Physical  characteristics  of  a missile  are: 

1.  Diameter  (in. ),  variable  name  D! 

2.  Weight  (lb. ),  variable  name  W 

3.  Max.  flipper  deflection  (deg. ),  variable  name  AG. 
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Other  Characteristic! 

Other  characteristics  are: 

1.  Lifespan  (sec. ),  variable  name  TMAX 

2.  Launch  velocity  (It. /•«. ).  variable  name  X(Z) 

J.  Max.  nipper  travel  <d«g. ).  variable  name  FMAX. 

Information  not  shown  but  required 

Information  not  shown  but  required  is: 

1.  Mle.Ue  HU  radlue  need  in  part  to  determine  probability  of  bit  of 
missile,  variable  name  RMK 

2.  Missile  name,  variable  name  MISTYPE 

3.  Blur  circle,  variable  name  BLURC 

4.  Minimum  detectable  Irradiance.  variable  name  HMIN. 

FLARE  FILE 

Tb.  flare  information  ha.  t»o  cat.,. Heat  4 

other  cbaracterUtice.  Tb.  flare  file  eubroutln.  1.  pr.e.nted  a.  Table  1-4. 


Table  1-4.  Flare  file  subroutine 
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i I 


: i 
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Not.  that  the  program  can  handle  pyrophoric.  U they  are  modelled  a.  .pedal 
flare.,  e.g.  .hort  rla.  time,  ahort  born  time,  high  peak  Int.n.lty  and  high 
drag  coefficient.  However,  the  program  may  require  aome  development  to 
more  fully  take  account  of  the  burning  and  aerodynamic  char.cterl.tlc.  of 
thl*  type  of  countermeasure. 
physical  Characteristics 

Physical  characteristics  are: 

1.  Diameter  (in. ),  variable  name  DF 

2.  Length  (in. ),  variable  name  XL 

3.  Weight  (lb. ),  variable  name  WFO 

4.  Grain  Weight  (lb. ),  variable  name  WG. 

Other  Characteristics 

Other  characteristics  are: 

1.  Burn  time  (sec. ),  variable  name  TBURN 

2.  Reference  intensity  (watts /ste rad),  variable  name  REF 

3.  Spectral  band  constant*,  variable  name  FIBAND 

4.  Drag  coefficient,  variable  name  CDF 

5.  Dispenser  ejection  velocity  (ft/sec),  variable  name  VFLARE 

6.  Flare  type  (name),  variable  name  IFTYPE 

7.  Type  of  surface  area,  variable  name  MK. 

TARGET  FILE 

Target  information  has  three  categories:  physical  characteristics. 
Initial  condition,  and  other  char.cterl.tlc..  Table  1-5  give,  the  target  file 
subroutine. 


Table  1-5.  Target  file  subroutine 

ivuouuw  r«/rt  o*»«i  *>i  k.*»»i»a  at/iim  i«.ia.^r 


fU'OritlttHWUli'i'OUMIM'Ull 

•r  k-  fr> 

m it  i*  a i*  i)  .*  u i*  <j  J J o j a it  tt  j ii 


Phytcal  Characteristics 


Physical  characteristic*  are: 

1.  Longitudinal  distance  from  tailpipe  to  tip  of  tail  (ft),  variable 
name  XB 

2.  Longitudinal  distance  from  tailpipe  to  nose  (ft.),  variable  name  XN 

3.  Wingspan  (ft. ),  variable  name  ZS. 

Initial  Condition 

Initial  condition  is: 

1.  Maximum  aircraft  turn  (g's),  variable  name  FMG 

2.  Maximum  aircraft  forward  acceleration  (g's),  variable  name  AM 

3.  Maximum  aircraft  speed  (mach),  variable  name  VM 

4.  Target  altitude  (ft. ),  variable  name  XP(8) 

5.  Target  velocity  (ft/sec.),  variable  name  X(7). 

Other  Characteristics 

Intensity  as  a function  of  polar  angle  are  tables 

1.  Polar  angle  (deg. ),  variable  name  PANG 

2.  Intensity  (kw/sterad),  variable  name  RINT 

3.  Target  type  (name),  variable  name  IAC. 

ATMOSPHERE  FILE 

The  last  file  to  be  discussed  is  the  atmospheric  spectral  transmittance 
tables.  These  tables  are  a function  of  range,  altitude,  temperature,  and 
band  region  for  target  and  flare. 

1.  Range  (ft. ),  variable  name  RNG 

2.  Atmospheric  spectral  transmittance  of  target,  variable  name  TAUT 

3.  Atmospheric  spectral  transmittance  of  flare,  variable  name  TAUF. 
Table  1-6  gives  the  atmospher  file  subroutine. 
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Table  1-6.  Atmosphere  file  subroutine 
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IDSTIUST  IT 

IDSTIUST  1* 
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IDSTIUST  fl 

IDSTIUST  tl 

IDSTIUST  tt 

IDSTIUST  IS 


DATA  SOURCES 


J A Ratkovic , K.  Nishimcto  'TRCM  Simulation  Study  ,Uj  ••  Quarteriy 
Progress  Report  prepared  by  Hughes  Aircraft  Company,  Culver  City, 
Calif,  for  Air  Force  Avionics  Laboratory,  Wright- Patterson  Air 
Force  Base,  Ohio,  July  1973 

IRCM  Simulation  Study  (U),  Feb.  1972,  Volume  H Hughes  Aircraft 
Co.,  Culver  City,  CA,  SDN  G-5812 

Fighter -Launched  Missiles  (Trends),  Eurasian  Communist 
Countries  (U),  Dec.  1971,  Defense  Intelligence  Agency,  Report 
No.  T65-09-26B.  SDN  J-58432 

DAWN  (Develop  Attack  Warning  Needs)  < U ) . FI nal  R e p ° rt.  Nov . 1972, 
General  Research  Corps.,  Santa  Barbara,  CA,  SDN  G-61305 

Foreign  Material  Exploitation  Report,  Interim  Report,  Project 
Graduation  Level  (U),  July  1972,  Missile  Intelligence  Agency,  U.S. 
Army  Missile  Command.  Redstone  Arsenal,  Alabama,  SDN  G-6135 

Shoulder  Fired  Surface  to  Air  Missile  System  Comparison 
Summary  (U).  June  1972.  Mlcom  Foreign  Intelligence  Office. r;  61358 
Army  Missile  Command,  Redstone  Arsenal,  Alabama,  SDN  G-61358 

AIM-9D  Simulation  Parameters,  Gene  Yowkl",  TechnUal  Note 
4055-2-68,  U.S.  Naval  Weapons  Center,  China  Lake,  Calif.,  Dec.  196 

Hughes  AIM-4D  Aerodynamic  Data,  SRS-585,  Revised  1 January  1965 

Assessment  of  Aerodynamic  Studies  of  Foreign  Tactical  Missiles, 
Leroy  Spearman  and  Charlie  Jackson,  Jr. , NASA  Langley  Research 
Center,  Hampton,  Va. , Feb.  1971 

ANAB  Missile  Wing  Evaluation,  FTD-CW-09-4-70,  Feb.  1970. 

Report  on  ASH  Alr-to-Air  Missile  Exploitation,  Ministry  of  Defense, 
TM  88169,  Dec.  1969. 

MAVWEPS  Report  TN  4063-233,  AIM-9L  Wind  Tunnel  Test  Report  (U), 
Oct.  19”2 

PMS  12AD44- 1/2174,  Final  Stability  »nd  Control  Report  for  the 
AIM-54A  Missile  (U),  27  April  1973 

Radiometric  Data  and  Mission  Profile,  Dept,  of  the  Air  Force, 
Headquarters  Aeronautical  Systems  Division,  Wrlght-Patterson  AFB, 
Dec.  1972,  SDN  G-61377 

B-52  Infrared  Radiation  Patterns,  Beoing  Co. , Wichita  Division, 

July  1968,  SDN  J-5V090 
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2.  PROGRAM  INITIALIZATION 


This  portion  of  the  simulation  program  initializes  all  program  constants, 
launch  geometry  variables,  and  determines  the  missile  level  angle  computations. 
See  Figure  2-i.  Section  2.1  describes  the  initialization  procedures,  and 
Section  2.  2 the  lead  angle  computations. 

2. 1 PROGRAM  INITIALIZATION 

In  the  missile,  target,  flare  simulation  program,  basic  parameters  have 
to  be  defined,  and  the  program  involved  in  doing  this  is  initial.  The  three 
categories  to  be  discussed  will  be  definitions  of  constants,  initialization,  and 
launch  geometry.  A block  diagram  of  these  computations  is  shown  in  Fig- 
ure 2-2.  See  Table  2-1  for  program  listing. 

Definitions  of  Constants 

The  following  constants  are  defined  in  the  simulation  program: 

1.  Gravity,  (ft/sec. ) variable  name  G,  32.  2 
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Table  2-1.  Initial  subroutine 
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Table  2-1  (Continued) 
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3.  Air  density  at  tea  level,  variable  name  RHOZ  {siuge/ft  ) and 
coefficient  of  exponential  variation  with  altitude,  variable  name 
CZ  (ft*1) 

4.  Atmospheric  density  aa  a function  of  altitude,  variable  name  RHO 

5.  Speed  of  sound  as  a function  of  altitude,  variable  name  FMACH 

6.  Surface  area,  variable  name  SFB  (ft^). 

Atmospheric  density  is  given  as: 

RHO  = RHOZ  ♦ EXP  (-CZ  * ALTITUDE) 

RHOZ  * STD  air  density  (sea  level) 


where, 

ALTITUDE  = XP(8)  = TARGET  ALTITUDE 

Speed  of  sound  is  a function  of  altitude.  Altitude  of  missile  and  target 
are  averaged  and  checked  whether. 


Altitude  > 36152,  FMACH  = 968.1 


IF 


* 


Altitude  s 36152,  FMACH  = 1116.0  - 0.0041  * Altitude 

• \ 

Initialization 

Pertinent  missile,  target  and  flare  parameters,  including  position, 
velocity,  and  acceleration  are  initialized.  Certain  parameters  previously 
defined  by  data  file  input  are  discussed  in  Section  1.  They  are: 

1.  Target  range  (ft.),  variable  name  X{6) 

2.  Target  velocity  (ft/sec.),  variable  name  X(7) 

3.  Target  altitude  (ft.),  variable  name  XP{8) 

4.  Missile  altitude  (ft.)*  variable  name  XP(3). 
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Other  required  initial  conditions  are: 

1.  Elevation  angle  as  a function  of  a random  number  (radians), 
variable  name  AV 

2.  Flare  time  between  salvos  (sec.),  variable  name  TF 

3.  Flare  dispenser  ejection  velocity  components  (ft.  /sec.),  variable 
name  VXF,  VYF,  VZF 

4.  Minimum  flare  velocity  (ft.  /sec.),  variable  name  VTERM. 

Minimum  flare  velocity  is  precomputed  here  to  be  used  later  as  a check 
on  radiant  intensity.  This  velocity  is  based  on  flare  burn  time,  surface  area, 
weight,  and  drag. 

Laun^.t  Geometry 

Options  available  for  missile  launch  conditions  are: 

1.  Input  azimuth  angle  (AH)  and  horizontal  range  (RXY)(in  which  the 
program  selects  the  elevation  angle  (AV)  on  a random  basis. 

2.  Input  azimuth  angle,  horizontal  range,  and  elevation  angle. 

3.  Input  azimuth  angle,  horizontal  range,  and  vertical  displacement 
(Vdisp). 

2.  2 LEAD  DETERMINATION 

This  subroutine  computes  initial  velocity  components  for  the  missile  at 
* launch.  Several  alternate  launch  modes  are  available.  Different  launch  modes 
may  be  used  in  the  elevation  and  azimuth  planes.  See  Figure  2-3  for  a block 
diagram  description. 
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Figure  2-3.  Lead  subroutine  block  diagram 

\ 

The  following  laws  are  available 

1.  Lead  collision  with  maximum  lead  limit. 

2.  Pursuit 

3.  Visual  lead  pursuit  with  maximum  lead  limit. 

4.  Pursuit  with  super-elevation  angle  (elevation  plane  only). 

Lead  collision  launch  is  based  on  attempting  to  put  the  missile  bn  a 
collision  course  for  a missile  velocity  of 

VL  + .'V, 

where  is  launch  velocity,  and  AV  is  an  incremental  velocity. 
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In  pursuit  launch,  the  missile  is  launched  on  a line  directly  toward  the 
target.  This  mods,  is  used  when  the  launcher  or  the  missile  does  not  have 
capability  for  lead  launch. 

In  visual  lead  pursuit,  a lead  angle  is  estimated  by  the  person  launching 
the  missile.  The  value  of  lead  is  generally  restricted  to  a small  angle  in  this 
mode.  The  lead  angle  nL  is  computed  from 

Sin  3 K sin  A^., 

where  is  target  aspect  angle,  and  K is  an  empirical  constant -0.  5. 

Often  the  missile  is  launched  in  a trajectory  above  the  target.  The 
angle  above  the  angle  of  launch  is  called  the  super-elevation  angle.  The 
superelevation  angle  is  input  as  a constant,  with  a limit  on  the  total  elevation 
angle  of  launch.  A program  listing  is  contained  in  Table  2-2. 

. V 
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Table  2-2.  Lead  subroutine 
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3.  DYNAMICS 

Figure  3-1  Is  a block  diagram  of  the  major  component*  of  the  simulation 
program  and  shows  the  dynamics  computations  to  be  performed.  Basically, 
the  dynamics  portion  of  the  simulation  program  determines  the  X-,  Y-, 
and  Z -components  of  acceleration,  velocity,  snd  position  for  for  the  missile, 
target,  and  flare(s). 

The  forces  controlling  the  missile  trajectory  are  thrust,  chord,  com- 
manded, and  gravity.  The  missile  velocity  and  position  are  computed  by 
Integrating  the  total  acceleration. 

The  target  Is  considered  to  fly  nominally  straight  and  level,  l.  e. , no 
maneuver.  However,  the  program  does  allow  for  three  maneuver  options: 

(1)  turn  In  any  direction,  (2)  straight  acceleration,  and  (3)  turn  with 
acceleration. 

When  the  flare(s)  Is  deployed,  the  flare  deployment  strategy  controls 
how  many  are  deployed,  In  what  direction  they  are  deployed,  and  how  often 
they  are  deployed.  The  forces  which  govern  the  flare  motion  are  drag  and 

gravity. 

3.1  MISSILE  DYNAMICS 

The  missile  equations  of  motion  are  governed  by  the  four  forces  shown 
In  Table  3-1.  This  table  also  lists  the  X-,  Y-,  and  Z-  components  of  each 
force.  Figure  3-2  shows  a block  diagram  of  the  computations  performed  In 
this  portion  of  the  program.  Basically,  the  missile  dynamics  are  updated 
as  follows: 

1.  The  almpolnt  position  Is  fed  Into  the  gyro  subroutine  from  which 
gyro  position  and  rate  are  output. 

2.  These  gyro  rates  and  positions  are  then  fed  Into  a subroutine  which 
simulates  the  guidance  unit  of  the  missile  and  computes  the  accelera 
tlon  components  which  are  to  be  commanded  by  the  missile. 

3.  The  acceleration  components  due  to  thrust  and  chord  forces  are 
subsequently  computed  and  resolved. 
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Table  3-1.  Force*  acting  on  the  miaaile 

THtUST 


A-  NAVtOATION  f.MAMITn  * m «• 

4.  The  mU.ile  acceleration  component,  due  to  thru.t.  chord, 
commanded  and  gravity  force,  are  .ummed  to  obtain  the 


and  Z-mi..ile  acceleration  component*. 

5.  Tb...  component*  nr.  then  lnt.gr.ted  to  obtein  the  mle.ll.  velocity 


component*. 

6.  Finally,  the.e  velocity  component,  are  .ub.equently  Integrated  to 
obtain  miaaile  position. 


The  remainder  of  thl.  Section  de.crlbe.  the  ,yro  po.ltlon  and  rat.  com- 

natation,  the  thru.t  computation.,  the  chord  force  computation.,  and  the 
patationa,  tn  „ ■>  _?  -how*  the  portion 

commanded  acceleration  computation.  In  detail.  Table  3 *.«  P 

of  the  main  program  vdtlch  Involve,  the  ml..U.  dynamic,  computation.. 
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Figure  3-1.  Dynamics 
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Table  3-2.  Executive  routine  ehowing  mlsslle/dyn&mlcs  computation* 


(Table  3-2,  continued) 
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(Table  3-2,  concluded) 
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Gyro  Position  and  Rate  Computations 

Because  there  ia  a time  lag  associated  with  the  gyro  and  the  forward 
tracking  loop,  the  aimpoint  location,  determined  in  aubroutine  aimpoint,  la 
not  tracked  preciaely  by  the  miaalle.  The  actual  value  of  4»  (aimpoint  rate) 
which  ia  requ’.  ec  to  command  proper  acceleration,  ia  not  fed  into  the  guid- 
ance unit  of  the  miaalle  until  aeveral  time  conatanta  later.  Thla  eftect  la 
modeled  in  the  program  by  a aimple  one  time  conatant  delay,  which  ia  given 

by: 

1 +oW  4 (^)  (+  ‘ *oU> 

where, 

A = integration  atep  aize 

e 

^new  * Precent  «yro  rite 

= previoua  gyro  rate 
* aimpoint  LOS  rate 

Tc  * time  constant  of  the  forward  tracking  loop 

Since  the  almulatton  uaes  a two-plane  geometry,  a computation  for 
ia  made  for  both  the  horizontal  and  vertical  plane  ualng  the  equation 
described  above.  If  the  gyro  rate  computed  from  thla  equation  becomes 
greater  than  maximum  precession  rate,  it  ia  set  to  the  maximum  value. 
The  gyro  position  is  obtained  by  an  integration  of  the  rate.  Figure  3-3 
shows  the  computations  performed  in  this  aubroutine. 

In  addition  to  computing  gyro  rate  and  position,  this  subroutine  also 
checks  the  miaalle  FOV  and  seeker  look  angle  limits.  If  either  of  these 
limits  are  exceeded,  the  program  goes  into  an  abort  mode  and  the  point  of 
closest  approach  ia  subsequently  computed  along  with  the  probability  of  hit. 
Table  3-3  contains  a listing  of  this  subroutine. 
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Figure  3-3.  Gyro  position  and 
rate  computations 
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Table  3-3.  Gyro  position  and  rate  computation  subroutine 
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Thrust 


The  propulsion  system  of  a missile  is  completely  defined  by  thrust- 
time history,  specific  impulse  (delivered)  and  motor  weight  drop  (if  applic- 
able) data.  Figure  3-4  shows  the  computations  performed  in  subroutine 
thrust  and  Table  3-4  contains  a listing  of  this  subroutine. 

The  thrust-time,  motor  weight  drop -time,  and  specific  impulse-time, 
profiles  are  stored  tables  in  the  program  which  are  read  in  as  part  of  the 
missile  flic  data.  The  values  of  these  variables  are  then  found  by  means  of 
a table  lookup. 

It  is  necessary  to  formulate  this  force  into  its  X-,  Y-,  and  Z-accelera- 
tion  components.  Table  3-5  shows  the  computational  procedure  for  performing 
this  operation  and  Figure  3-4  shows  the  equations  used  in  implementing  this 
component  resolution,  with  the  (G /W  nflW)  factor  accounting  for  the  conversion 
of  force  to  acceleration. 

The  thrust  subroutine  in  addition  to  computing  the  thrust  components 
also  computes  the  missile  weight.  The  change  in  missile  weight  during  the 
thrust  period  is  given  by: 

W = W . . - (thrust/S)  A - W 
new  old  o 


where 

W 


new 


W 


old 

A 


W 


new  missile  weight 
old  missile  weight 
integration  step  size 
specific  impulse 

motor  weight  drop  at  end  of  boost  period  (if  applicable) 
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TIME  1. 


Thrust  computations 


Table  3 -4.  Thruat  aubroutine 
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Chord 


The  magnitude  of  chord  force  is  given  by  the  expression 

l/2pCcn(DM/2)2V2 


where  i 

p s atmospheric  density 

C = chord  force  coefficient  (a  function  of  mach  number) 
c 

tt(D.-/2)2  = missile  cross-sectional  area 
M 

V s missile  velocity 

The  atmospheric  density  (p)  is  modeled  as  an  exponential  function  of 
altitude  and  computed  only  once  in  subroutine  initial. 

The  chord  force  coefficient  Is  a function  of  mach  number  and  is  found 
in  the  program  by  me*ns  of  a table  lookup. 

The  missile  velocity  needed  to  calculate  the  chord  force  is  computed  in 
the  range  and  range  rate  computation  subroutine. 

The  chord  force  is  directed  opposite  the  thrust  vector  and  therefore 
has  the  same  X-,  Y-,  and  Z-unlt  vector  components  as  the  thrust  vector. 

Figure  3-5  shows  the  overall  computations  of  subroutine  thrust  and 
Table  3-6  contains  a listing  of  this  subroutine. 
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Chord  acceleration  computation 


Table  3-6.  Chord  eubroutlne 
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Commanded  Acceleration 

Commanded  Force 

The  commanded  force  results  from  the  horizontal  and  vertical  guidance 
commands  generated  in  the  missile.  The  commanded  force  is  based  on  a 
proportional  navigation  system  with  navigation  parameter.  It  is  assumed 
that  the  missile  has  its  control  surfaces  or  flippers  biased  to  account  for 
any  lift  force  and  ti  fly  straight  and  level  during  normal  flight  conditions. 
The  commanded  acceleration  in  the  vertical  and  horizontal  plane  is  given 
by  the  expression: 


(horizontal) 


(vertical) 


navigation  parameter 

gyro  rates  - horizontal  and  vertical 

gyro  position  - horizontal  and  vertical 

angle  of  attack  . horizontal  and  vertical 

missile  body  angles  - horizontal  and  vertical 

aimpoint  range  rales  - horizontal  and  vertical 

gravity  bias  term 


The  derivation  of  this  guidance  law  is  given  in  the  following  subsection. 

To  account  for  gravity  bias  missile  systems,  i.e.  , missiles  which 
have  their  horizontal  control  surfaces  biased  in  such  a manner  as  to  effec- 
tively null  out  gravity,  the  vertical  commanded  acceleration  includes  a 
gravity  bias  term,  G*BIAS. 
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BIAS  It  the  variable  which  control!  the  amount  of  gravity  bias  the 
mittile  it  to  have. 

To  compute  the  commanded  acceleration,  it  la  neceaaary  to  determine 
• • • 

r^y  and  r^  *n  terma  of  r (the  elating  rate  along  the  LOS).  Table  3-7 
thowt  thia  computation. 

Figure  3-6  thowt  the  commanded  acceleration  computation  and  also 
ahowa  that  commanded  acceleration  ia  aerodynamlcally  and  structurally 
limited. 


The  commanded  force  it  limited  aerodynamically  to  be  lett  than  the- 
maxlmum  lift  force  given  in  the  same  table  to  be  l /2*p*Cj^^2£*ir*(Dj^/2)^<‘V^ 
CnMAX  **  * ^unct^0n  °*  rnach  number  aa  incidated  in  the  figure. 

A limit  is  set  in  the  missile's  autopilot  or  guidance  unit  to  prevent  over- 
maneuvering  against  the  target.  This  ia  a g-llmit  and  la  designed  by  the 
term  A5  in  tho  program. 

The  g-limit  it  a limit  internal  to  the  missile  whereat  the  aerodynamic 
limit  itan  external  limit.  At  any  given  time  and  for  any  given  missile,  only 
one  or  the  other  constraint  will  be  dominant.  These  two  limits  represent 
constraints  on  commanded  acceleration  X(5),  XP(5))  and  are  also  the  only 
state  constraints  in  the  program. 

The  missile  does  not  respond  instantaneously  to  the  commanded  force. 
There  is  a delay  associated  with  time  for  target  information  to  go  through 
the  signal  processor  and  guidance  unit  and  finally  to  reach  the  control  surface 
actuators.  This  delay  is  also  modeled  as  a one  time  constant  delay.  The 
resulting  equations  are: 


X(5), 


new 


= X(5) 


old 


(G 1 - X(5) 


old) 


XP(5) 


new 


- XP‘S>oW  <01'  - XP<S»old> 
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Tab!  e 3-7.  Range  rate  vector  component* 


From  geometry, 


*x  = *xy  * co“  M “ \z  * COS 


rxz  = rxy  * cos  (4»)/cos  (41' ) 


Now, 


•2  *2,  • 2 , >2 
r a r + r + r 
x y z 


rxy2  * cos2  (41)  + rxy2  * sin2  (4)  + r^2  * sin2  (+') 


• 2 . . 2 . 2 , 
r + r + sin  (tL1) 
xy  xz  ' 


rxy2  + r xy2  * sin2  (+•)  * cos2  (^)/cos2  (41) 


xy 


1 + cos2  (4)  * tan2  (41) 


|T72 


Figure  3-6.  Commanded  acceleration 
. computations 


1 


where, 

G1  = commanded  acceleration  required  (Horizontal  plane) 

I 

G1  = commanded  acceleration  required  (vertical  plane) 

X(5)  = actual  acceleration  at  the  control  surfaces 

XP(5)  = actual  acceleration  at  the  control  surfaces  (vertical) 

T = missile  time  constant 

s 

The  commanded  fo~ce  is  in  a direction  normal  to  the  thrust  vector. 

It's  X-,  Y-,  and  Z-vector  components  are: 

Ax  * - sin  (V  + a)  * X(5)new  - sin  (y  + <* ) x XP(5)ftew 
Ay  = cos  (Y+«)  * X(5)new 

Az  = cos  (V  + «)  * XP(5)riew 

Table  3-8containsa  listing  of  this  subroutine. 

The  missile  angle  of  attack  is  also  computed  in  this  subroutine.  Oue 
angle  of  attack  (alpha  angle)  i3  computed  for  the  missile  in  the  vertical  plane 
and  one  in  thehorlzontal  plane.  Curves  of  alpha  angle  as  a function  of  Mach 
No.  and  normal  force  coefficient  are  generally  available  from  missile  speci- 
fications. This  data  is  tabularized  into  a two  dimensional  array  of  alpha 
as  a function  of  mach  number  and  normal  force  coefficient  for  use  in  the 
program.  The  procedure  for  determining  alpha  in  each  plane  is  to  compute 
the  normal  force  coefficient,  based  on  the  achieved  accelerations,  and  the 
mach  number  and  then  do  a table  lookup  for  alpha. 

Guidance  Law 


In  summary,  the  guidance  law  implemented  in  infrared  missile  auto- 
pilots is  generally  of  the  form 


Table  3-8 


Commanded  acceleration  subroutine 
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where. 


A^,  * commanded  acceleration  normal  to  body  axis 
4>  2 juOS  rate  • 

Kg  = navigation  gain 

Because  the  navigation  gain  (Kg)  varies  considerably  with  tactical  con- 
ditions (altitude,  launch  speed,  target  speed,  etc.)  it  is  difficult  to  decide 
upon  a fixed  value  for  this  constant  when  working  with  one  missile  only. 

When  many  missiles  are  to  be  evaluated,  as  is  the  case  of  the  missile, 
target,  and  flare  simulail  n,  this  guidance  law  becomes  too  difficult  to  imple- 
ment and  can  lead  to  a .*e  variation  in  missile  trajectory.  In  order  to 
avoid  this  problem,  the  missile,  target,  and  flare  simulation  uses  an  ideal 
guidance  law  given  by 

Ari 

AT  * coa  (y  - 4*  + ® ) 

where  these  parameters  are  defined  in  Table  1.  It  should  be  pointea  out 
that  this  law  is  impossible  to  implement  in  IR  missile  hardware  due  to  the 
fact  that  missile  and  target  range  rate  information  is  required.  However, 
it  does  enable  the  simulation  to  remain  generic  and  to  evaluate  missile  per- 
formance under  ideal  conditions. 

The  guidance  law  of  the  form  Kg4<  can  also  be  implemented  in  the  simu- 
lation, particularly  when  statistical  variations  in  missile  trajectory  are 
desired.  This  can  be  accomplished  easily  by  varying  the  navigation  gain 
over  its  range  of  values. 

The  commanded  missile  acceleration  normal  to  the  iine-of-sight  (LOS) 
for  ideal  proportional  navigation  is  given  by 

An  = AVC  + 

Figure  3-7  shows  the  reference  geometry  and  Table  3-9  defines  the 
glossary  of  terms. 


Figure  3-7.  Missile  and  target 
reference 
geometry 

l.i  actual  practice,  the  commanded  acceleration  is  directed  normally  to  the 
body  axis  and  given  by 


, Ar»ji 

AT  “ COS  Of  cos  (®  - 


This  relationship  is  the  equation  for  Implementing  the  proportional  naviga- 
tion law  in  the  ideal  case,  and  is  presently  the  equation  used  in  the  missile, 
target,  and  flare  simulation. 

From  the  infrared  missile  hardware  point  of  view,  this  equation  is 
impossible  to  implement  because  range  rate  information  is  not  available  to 
thr.  missile.  The  equation  generally  implemented  in  the  autopilot  of  infrared 
missiles  is 

At  = KJ  (3-8) 

T g 

where 

AV 

Kg  ■ c..(vA—  <3-” 
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Tabic  3-9.  Glossary  of  terms 

V , V » inertial  reference  axes 
x y 

Vj  * target  velocity 

= angle  between  target  velocity  and  inertfral  axis 

V b missile  velocity 
m 1 

y s angle  between  missile  velocity  and  inertial  axis 

4*  = anglr  between  LOS  and  Inertial  axis 
• 

s LOS  rate 
<r  s angle  of  attack 
T s missile  body  axis  direction 

e 

v s range  rate  along  LOS 
Vc  s missile /target  closing  velocity 

3 missile  acceleration  normal  to  LOS 
A .j.  s missile  acceleration  normal  to  body  axis 

A r navigation  parameter 
K 3 navigation  gain 

I 2 

The  parameter  is  set  before  launch  and  is  a function  of  altitude,  target 
speed,  and  launch  speed.  This  constant  is  normally  chosen  to  be  three  times 
the  maximum  closing  velocity  achieved  by  the  missile  during  flight.  This 
ensures  that  the  minimum  value  of  the  navigation  parameter  will  be  3,  The 
values  during  flight  will  be  higher,  ranging  up  to  5 or  more  at  the  end  of 
flight.  The  constant  value  of  4 used  in  the  simulation  is  an  approximation 
of  the  average  value  during  an  actual  flight.  This  is  somewhat  optimistic 
from  the  missile  point  of  view. 

3.  2 TARGET  EQUATIONS  j 

< 

The  program  allows  for  four  different  types  of  flight  paths  to  be  flown 
by  the  target  aircraft.  These  include:  straight  and  level,  circular  turn, 
straight  acceleration,  and  turn  and  tangential  acceleration. 

All  target  motion  is  considered  to  take  place  only  in  the  horizontal  X-Y 
plane,  and  the  equations  governing  each  Right  path  are  given  in  Tables  3-10 
through  3-13. 
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Table  3-10.  Straight  and  level  flight 


Equations 

V 

V 

ii 

K 

H 

> 

VT 

co»  yt 

VT 

*y 

VT  Y = 

VT 

• In  yt 

) 

VTz  = 

0 

Table  3-11.  Circular  turn 
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Table  3-13.  Turn  and  acceleration 


The  target  altitude  and  velocity  ia  an  initial  input  into  the  program.  The 
target  heading  relative  to  the  inertial  axes  la  determined  in  the  program  by 
miaaile  aapect  angle  at  launch.  The  additional  aircraft  data  required  to  effect 
the  tnaneuvera  includea  the  maximum  number  of  g'a  the  aircraft  can  pull  in  a 
turn*  the  maximum  aircraft  acceleration  forward,  and  the  maximum  apeed  it 
can  obtain. 

Figure  3-8  ahowa  the  baaic  equationa  uaed  In  the  program  to  simulate 
target  motion.  Target  maneuvers  are  controlled  through  NAMELIST  input 
variables  TM,  AM,  FMG.  The  time  of  initiation  of  maneuver  ia  set  by  TM. 
Linear  acceleration  AM,  and  target  turns  FMG,  are  input  as  the  magnitude 
of  the  acceleratlon(s)  in  terms  of  the  number  of  g'a.  A positive  sign  for 
FMG  will  generate  a right  turn  and  a negative  sign  a left  turn.  Table  3-14 
contains  a listing  of  this  subroutine. 
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Table  3-14 


Target  motion  subroutine 
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3.  3 FLARE  CONTROL  AND  DYNAMICS 

Figure  3-9  shows  a block  diagram  of  the  major  computations  performed 
in  this  section.  Basically,  the  flare  control  portion  of  this  problem  is 
responsible  for  determining  the  flare  status  (not  dispensed , burning,  or 
extinguished),  the  times  at  which  flares  are  ejected,  and  the  initial  velocity 
and  position  of  the  flares  at  ejection.  The  flare  dynamics  updates  the  flare 
trajectory  for  those  flares  which  have  been  dispensed  and  are  still  burning. 

Flare  Control 

The  flare  status  as  indicated  in  Figure  3-10  is  controlled  by  the  variable 
NM(k),  with  NM(k)  = l indicating  that  the  flare  has  not  yet  been  dispensed, 
NM(k)  = 2 indicating  that  the  flare  has  been  dispensed  and  is  still  burning, 
and  NM(k)  = 3 indicating  that  the  flare  has  been  extinguished.  The  flare 
initial  deployment  strategy  is  also  under  the  control  of  this  portion  of  the 
program.  Figure  3-10  shows  that  flares  can  bi  dispensed  as  a function  of 
(l)  time  after  missile  launch,  (2)  missile  and  target  range,  *3)  time  to  go, 

and  (4)  missile  and  target  range  rate. 

Once  a flare  has  been  deployed,  the  flare's  Initial  velocity  is  computed 
based  on  the  target  velocity  plus  the  dispenser  ejection  velocity.  This 
includes  the  ejection  direction.  The  flare's  position  is  computed  based  on 
target  position  plus  dispenser  location  relative  to  engine(s)  position.  A 
listing  of  this  subroutine  is  contained  in  Table  3-15. 


FLAM  CONTROL  PU«  DYNAMICS 


1.  DETERMINE  EACH  FLAM  STATUS 

• NOT  DISPENSED 

• EJECTED  ANO  MINING 

PCS  THOSE  PIARES  EJECTED  AND  STIll  MINING 

PLARE 

• EXTINGUISHED 

acceleration 

1.  DETERMINE  PLARE  .‘SPACE  AREA 

VELOCITY 

MISSlUAAKSn 

».  DETBMINE  WHEN  TO  DEPLOY  EACH  FLARE 

i.  DETERMINE  PLARE  WEIGHT 

POSITION  _ 

• TIME  AFTER  LAUNCH 

■* 

• RANGE 

1.  COMPUTE  DRAG 

• TIME  TO  GO 

4.  DETERMINE  VELOCITY  AND  POSITION 

• RANGE  rate 

S.  DETERMINE  INITIAl  POSITION  AND 

VELOCITY  OP  EACH  FLARE 

DOIf-l,  N 

N • TOTAl  NUMIU  OF  FIXES 


TIME  (T) 

MISSILE/TAAGET  EANCE  id 
TIME  TO  CO  (TTG) 


IF. 
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SET  NM(ig  ■ J,  TfID  • Ti 

OPTION  2 

now  .no 

SETNMIlO-*,  yiQMi 

OPTION  3 
1 

1 CONSTANTS 

TF(K) 

• TIME  AFTU  LAUNCH  Vi  MEN 

FLAW  EJECTED 

•F00 

- UNGE  AT  WHICH  EJECTION 

OCCUES 

TTO(K) 

•TIME  TO  CO  WHEN  FIXE 

EJECTED 

FLAW  INITIAL  VELOCITY  COMKJTATION 


V^INITIAl'  *T  f11  ’ TF,K,j  *VDX 

V^imtial  - fT  - tf,k>]  *vdy 
V^INITIAL  ’ *T  [T*Tf"°]  *VDZ 


Table  3-15.  Flare  control  aubroutine 
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Flare  Dynamic! 

In  any  mlaaile,  target,  and  flare  almulatlon,  It  ia  vital  that  the 
trajectory  time  histories  of  the  flares  under  consideration  be  known. 
Tablo  3-16  shows  the  two  forces  acting  on  the  flare,  drag  and  gravity, 
as  well  as  the  magnitude  of  each.  The  following  equation  governing  the 
trajectory  time  history  can  then  be  written: 

VF  • -(7  PA.CD(0/WF)VFVF  ♦ o) 


where 

p s atmospheric  density 
Cp  « drag  coefficient 

* cross-sectional  area 
WF  * weight  of  flare 
Vp  = velocity  of  flare 
G = gravitational  constant 

The  major  mfflcultles  In  solving  this  equation  are  that  both  the  cross- 
sectional  area  and  the  weight  arc  time  dependent  functions,  while  the 
cross-sectional  area  Is  spatially  dependent  If  the:  flare  Is  tumbling  (fs  It 
normally  does). 

To  deal  with  the  spatial  orientation  problem,  It  Is  assumed  that  a cross 
sectional  area  averaged  over  all  surfaces  will  account  for  this  orientation 
of  the  flare  as  It  tumbles.  Then,  Ihe  cross-sectional  area  (SFB)  presented 
to  the  wind  stream  will  be  the  average  of  longitudinal  (A2)  and  axial  cross 
sectional  area  or 


A1  + A2 

— z — 


SFB 


Table  3-16.  Forces  acting  on  the  Kth  flare 


MAG 
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WF(IC)  - WEIGHT  OF  K*  FIA*E  5,(10  ■ SUIT  ACE  A1EA  OF  M«  FI  * * NOFMAL  TO  '.tlOCITY  VECTOt 

VF(W  • VnOCITY  OF  KH<  FUJI  C0  « MAO  COEFFICIENT 


OtAVITV 


Now  consider  cylindrical  flares  such  as  the  MK-46,  MK  -49,  and  ALA- 17. 
Assume  that  the  linear  burn  rate  is  constant,  and  that  the  area  of  the  cylinder 
walls  is  much  greater  than  that  of  the  ends  of  the  cylinder.  Thir  is  a good 
assumption  in  /lew  of  the  shapes  of  the  flares  being  investigated.  However, 
if  the  area  of  the  ends  of  the  cylinders  are  neglected  ancl  the  definitions  in  the 
sketch  below  are  used,  the  perimeter  surface  area  can  be  expressed  as 

SURFACE  AREA 


AMI  • 2 rfltl  /Ml 

On  the  assumption  that  the  linear  turn  rata  is  constant  and  that  the 
radius  of  the  cylinder  is  much  bhorter  than  .he  length  of  .he  cylinder,  it 
follows  that  dr/dt  = constant.  The  solution  .o  this  equation  is 

r(t)  = rQ(l  - t/tfl) 

\ 
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where 


r = Initial  radius  of  flare 

o 

t_  = burn  time 

D 

If  a similar  expression  for  f(t)  is  assumed,  which  is  not  necessarily 
the  exact  solution  for  this  term  but  one  that  should  result  in  an  adequate 
approximation,  the  expression  for  the  f(t)  becomes 

|(t)  = !0  (1  - t/tfi) 


For  cylindrical  type  flares, 

A l (t)  = rrr2(t) 


and 

A2(t)  = 2r(t)l(t)  w 

The  average  cross-sectional  area  is 

SFB(t)  . A1|t|  l * 2r(°  lltLj 

The  weight  of  the  flare  is  generally  specified  in  terms  of  both  a total 
weight  (WFO)  and  a grain  weight  (WG).  The  grain  weight  is  changing  as  a 
function  of  time.  If  it  i*  assumed  for  simplicity  that  the  grain  weight  changes 
with  a similar  expression  as  the  cross-sectional  area,  then  the  expression  for 
the  flare  weight  is  as  shown  in  Table  3-17.  This  tabic  shows  the  basic  com- 
putations the  simulation  uses  to  determine  the  flare  trajectory.  Figure  3-11 
shows  a flow  diagram  of  the  computations  of  this  subroutine  including  com- 
ponent resolution  of  all  acceleration,  velocity,  and  position  components. 

Table  3-18  con.....s  a listing  of  this  subroutine. 
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Table  3-17.  Drag 


Drag  acceleration*:  1/2  p Cp*"*Sp(k)V  /Wy(k) 

1.  2. 

1.  A1  * vrQ 

A2  = 2r  I w 
o o 

SFB  * (A1  + A2)/2 
SF(k)  = SFB  (1  - (T  - TF(k)/TBURN)2 

2.  WFtk)  = WFO  - WG(l  - (1  - (T  - Tp(k))/TBURN)2) 


3 


SUMAC!  AKA  computation 


/ 
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Table  3-18.  Flare  dynamic!  aubroutlne 
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4.  ANGLE  AND  RANGE  COMPUTATIONS 


In  this  acction,  the  missile  and  target  and  missile  and  flare  computations 
are  performed  to  determine  angle,  angle  rate,  range,  and  range  rate.  Sec- 
tion 4.  1 discusses  the  angle  computations,  and  Section  4.  2 discusses  the 
range  computations. 

4.  1 ANGLE  AND  ANGLE  RATE  COMPUTATIONS 

Tables  4-1  and  4-2  show  how  the  angles  between  the  missile  and  target 
and  their  corresponding  rates  shown  in  the  geometry  are  computed  along  with 
the  variable  name  associated  with  each  in  the  computer  program.  Subroutine 
MTGTANG  in  the  program  is  responsible  for  performing  these  computa- 
tions. See  Table  4-3  for  a listing  of  this  subroutine.  Figures  4-1  and  4-2 
show  the  computational  procedures  for  these  calculations  in  block  diagram 
form. 

Figures  4-3  and  4-4  describe  the  angles  and  angular  computations  for 
the  angles  between  the  missile  and  flare  flare).  Subroutine  MFLANG 
is  used  to  compute  these  angular  variables.  See  Table  4-4  for  a program 
listing. 


Angle  and  angle  rate  computation*  between  target  and  missile,  vertical  plane 


ZP(1 1)  [XP(6)*XP(9)-XP(7)*XP(8)1/[XP(6)*XP(6 

ZP(12)  ALPHAP+ZP(5)-ZP(3) 


Tabic  4-3.  Micelle  end  target  angle  and  angle  rate 
computations  cubroutinec 
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Missile  and  flare  angle  computations,  horizontal  plane 


Table  4-4.  Mi**lle  and  flare  angle  and  angle  rate 
computation*  subroutine 
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4.  2 RANGE  AND  RANGE  RATE  COMPUTATIONS 

Figure  4-5  show*  the  computational  procedure  for  determining  the 
mifcile  velocity  missile  acceleration,  miaalle  mach  number,  the  target 
velocity,  target  acceleration,  the  component  of  relative  position  velocity 
and  acceleration  between  the  missile  and  the  target,  and  the  relative  range 
and  range  rate.  In  addition  to  these  computations,  this  subroutine  also  cal- 
culates the  miss  distance  and  time-to-go  parameters.  There  is  a check  in 
the  computations  to  determine  if  the  range  rate  is  po stive  after  missile  thrust- 
ing is  terminated.  This  thrusting  termination  period  is  nominally  set  at 
5 seconds.  If  this  constraint  is  violated,  the  program  goes  into  an  abort. 
Tables  4-5  and  4-6  contain  listings  of  these  subroutines. 

Uigurs  4 ihow,  how  the  ran*e  between  the  missile  and  *n  arbitrary 
flare  (K  flare)  as  well  as  their  corresponding  rates  are  computed. 


/ 

I 


Table  1-5.  Range  and  range  rate  between 
mlcBilvJ  and  target  aibrcutines 
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Figure  4-6.  Mieeile  and  flare  range  and  velocity  computation* 


5.  TARGET  AND  FLARE  IR RADIANCE  COMPUTATIONS 


This  portion  of  the  program  determines  the  irradiance  from  the  target 
and  flare  at  the  missile  dome.  Section  S.  1 describes  the  target  irradiance 
computations,  and  Section  5.  Z describes  the  flare  irradiance  computations. 

5. 1 TARGET  IRRAPL'vNCE  COMPUTATIONS 

In  determining  the  irradiance  at  the  missile  from  the  target  there  are 
several  major  factors  which  must  be  determined  as  indicated  in  Figure  5-1. 

First,  the  aspect  angle  at  which  the  missile  views  the  aircraft  must  be 
determined.  This  is  found  from  the  dot  product  relationship. 

VT  • r = VT  r cos  6 
or 

XT  rx  + Yt  ry  + Z Trz  = VT  r cos  6 

where, 

V.p  = target  velocity 

r = missile/target  range 

X.J.,  Yj,  tj  = X,  Y,  Z components  of  target  velocity 
rX*  rY*  rZ  = Z components  of  missile/target  range 

0 = aspect  angle. 

The  target  intensity  data  is  stored  in  the  program  as  a function  of  polar 
angle.  Once  the  aspect  viewing  angle  has  been  determined  the  value  of  the 
target  intensity  is  found  by  means  of  a table  look  up  on  this  data. 

Data  on  atmospheric  attenuation  is  stored  in  the  program  as  a function 
of  range  and  black  body  temperature.  Based  on  the  aircraft  tailpipe  tempera- 
ture and  the  missile/target  range,  the  atmospheric  attenuation  is  found  by 
means  of  a table  look  up. 
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Figure  5-1.  Target  irradiance  computation 


/ 


The  lrradiance  from  the  target  it  then  computed,  at  Indicated  in 
Figure  5-1,  baaed  on  the  target  intensity,  atmospheric  attenuation  and 
missile /target  range.  A listing  of  this  subroutine  is  contained  in  Table  5-1. 
Notes  The  EPICS  program  consists  of  ASDIR  II  in  conjunction  with  the 
SPKINT  subroutine  and  M/T/CM.  As  we  have  just  seen,  M/T/CM  contains 
an  atmospheric  transmission  file.  Thus,  when  using  ASOIR  II  (or  auy 
other  program  which  generates  spectral  radiant  intensity)  it  is  Important  to 
use  a true  (unattenuated)  spectral  radiant  intensity. 

If  it  is  desired  to  use  a different  atmospheric  transmission  model,  this 
latter  model  must  be  used  (in  conjunction  with  target  temperature,  altitude, 
range  and  optical  waveband)  to  generate  a new  atmospheric  file  in  the 
M/T/CM  program. 


Table  5-1.  Target  lrradiance  subroutine 
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5.2  FLARE  IK  RADIANCE  COMPUTATIONS 

The  program  allows  for  two  options  in  determining  the  flare  Intensity 
profile  as  indicated  in  Figure  5-2.  Option  1 uses  a table  look  up  procedure 
for  determlng  the  flare  intensity  as  a function  of  time.  This  procedure  is 
generally  used  when  dynamic,  in  flight  Intensity  profiles  are  available  for 
altitude  and  wlndstream  conditions  at  or  near  those  being  considered.  When 
no  dynamic,  in  flight  data  exists,  option  2 can  be  used  to  generate  an  intensity 
profile  which  accounts  for  altitude  and  wlndstream  degradation  factors. 
Basically  this  model  consists  of  four  factors.  First,  a peak  flare  intensity 
value  (REFI)  measured  statically  and  referenced  to  the  1.7-2.  7p  band  is 
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Irr&dlance  computation* 


required.  Second,  the  variation  in  the  flare  ia  burning  surface  area  aa  a 
function  of  time  ia  required.  The  expression  shown  is  for  a cylindrical  type 
whose  linear  burn  rate  is  assumed  to  be  constant.  (See  reference  1. ) Next, 
the  factor 


exp 


V term  - Vp  (K) 
C 


accounts  for 


altitude  and  windstream  degradation.  This  is  an  empirical  expression  found 
to  curve  fit  measured,  dynamic  flare  data  quite  well.  The  term  V term 
represents  the  minimum  flare  velocity  and  is  determined  from  trajectoiy 
parameters  in  the  initialization  portion  of  the  program.  (See  also  Fifth 
Quarterly  Report,  IRCM  Simulation  Study).  Finally,  the  term  FIB  AND  is 
used  to  proportion  out  the  flare  energy  in  the  missile  band  being  considered 
relative  to  the  reference  band  (1.7-2.  7p). 

Once  the  flare  intensity  has  been  determined,  the  atmospheric  attenuation 
is  found  by  means  of  a table  look  up,  and  the  Irradiance  value  is  computed.  A 
listing  of  this  subroutine  is  contained  in  Table  5-2. 


Table  5-2.  Flare  irradiance  subroutine 
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6.  . AIMPOINT  DETERMINATION 


The  primary  functions  of  this  portion  of  the  program  are  to 

1.  Determine  which  infrared  sources  are  within  the  missile  FOV. 

2.  Determine  which  of  these  sources  within  the  missile  FOV  have 
irradiance  levels  above  the  minimum  detectable. 

3.  Determine  the  missile  aimf>olnt,  on  the  basis  of  the  missile's 
signal  processing  and  source  irradiance  levels  --  within  FOV  and 
above  minimum  detectable. 

Figure  6-1  describes  the  computational  procedure  for  determining  the 
aimpoint  in  block  diagram  form.  A further  detailed  description  of  these 
computations  is  contained  in  the  following  paragraphs.  A listing  of  this  sub- 
routine is  contained  in  Table  6-1. 

Table  6-2  shows  the  equations  used  to  check  which  ignited  flares  are 
within  the  FOV  of  the  missile  and  they  are  stored  in  an  array  Np(J)  for 
further  aimpoint  processing.  Similarity,  the  target  is  checked  to  see  if  it  is 
within  the  FOV  and  the  information  on  whether  it  is  or  not  is  stored  in  the 
variable  KT'. 

The  total  number  of  flares  in  the  FOV  is  indicated  by  the  variable  LT'. 
The  sum  of  the  variables  (LT'  + K’™)  Indicates  the  total  number  of  IR  sources 
in  the  FOV.  If  there  are  none,  then  the  program  will  go  into  an  abort  mode 
due  to  the  fact  that  there  are  no  Infrared  sources  in  the  missile  FOV. 

If  there  is  at  least  one  infrared  source  in  the  FOV,  the  program  deter- 
mine which  IR  sources  have  irradiance  levels  above  the  minimum  detectable 
by  the  missile.  The  flares  which  meet  this  criteria  are  stored  in  an  array 
N-,(J)  for  further  aimpoint  processing  with  their  total  number  ir.  the  array 

I 

being  indicated  by  the  variable  LT.  Similarily,  the  target  is  checked  to  see 
if  it  is  above  threshold  level  and  the  information  on  whether  it  is  or  not  is 
stored  in  the  variable  KT.  If  the  sum  of  these  variable  (LT  + KT)  equals 
zero,  then  the  program  will  abort  due  to  the  fact  that  there  are  no  IR  sources 
within  the  missile  FOV  above  threshold  value.  If  there  is  at  one  source  which 
meet*  this  criteria  the  program  will  then  determine  the  aimpoint. 
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Table  6-2.  Missile  FOV 
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The  aimpoint  type  can  be  baaed  on 

1.  geometric  centroid 

2.  ir radiance  centroid 

3.  maximum  irradlance 

of  IR  source,  within  the  FOV  above  the  minimum  detectable  irradiance  level. 
In  general,  con-scan.  FM  signal  processing  missiles  are  max  irradiance 
tracker,  and  spin-scan  AM  signal  processing  missiles  are  irradiance 

centroid  trackers.  . . , 

The  aimpoint  determine,  the  angle.  4»A.  »t'A.  angle  rate.,  +A.  A an 
the  rang.  rat.  r.  which  r.pre.ent  the  direction,  direction  rate,  and  the 
range  rat.  irom'the  ml. .lie  to  an  apparent  target  within  the  FOV.  The  equa- 
tion. u.ed  to  calculate  +A,  *A.  +X  *”d  'k  each  aimP°‘"*  C°ndi“°n 

are  shown  in  Table  6-3. 

The.,  aimpoint  variable,  arc  fed  hack  Into  the  dynamic,  portion  o I the 
program  to  determine  gyro  po.itlon  and  rate  and  ultimately  to  determine  the 

misoile  guidance  commands. 
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7.  CLOSEST  APPROACH  COMPUTATION 


Thl.  tub  routine  1.  celled  when  the  ml..lle  no  longer  ha.  a .ource  within 
the  aeeker  FOV.  The  ml..lle  and  target  trajectorle.  are  projected  foreward 
In  time,  a.tumlng  ml..lle  and  target  acceleration,  remain  con.tant,  at  the 
value,  before  lo..  of  tracking* 

Figure  7-1  *how.  a block  diagram  of  the  computation.  The  trajectory 
1.  prelected  (or weed  until  the  mle.ll.  pe.ee.  the  terg.t  (norm.l  t.rmlnetien). 
begin,  to  diverge,  hit.  the  ground,  or  exceed,  the  mexlmum  ml..U.  llletime. 

The  computed  time  to  go  (TTG)  1.  the  time  remaining  to  clo.e.t  approach 
aeeumlng  con.tant  ml..il.  and  target  velocltl...  It  1.  equal  to  the  proj.ct.on 
of  the  rang.  In  the  relative  velocity  direction,  divided  by  the  magnitude  of  the 
relative  velocity.  Table  7-1  give,  the  clo.e.t  approach  computation 
•ubroutlne. 
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Table  7-1.  Closest  approach  computation 
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8.  PROBABILITY  OF  HIT 


It  i.  necessary  when  evaluating  thousands  of  computer  runs  (1)  to  use 
the  probability  of  hit  (PH)  as  the  only  measure  of  effectiveness  in  the  simula- 
tion and  (2)  to  have  a simple  means  of  computing  PH  «o  as  to  keep  overall 
program  complexity  and  computation  time  to  a minimum.  The  approach 
taken  here  to  calculate  PH  utilises  the  following  assumptions  and  definitions: 

1.  The  missile  aimpoint  Is  located  at  the  geometric  centroid  of  all  the 
aircraft's  tailpipes  and  thus  the  point  of  missile  closest  approach 
to  the  aircraft  is  relative  to  the  tailpipe. 

2.  The  tailpipes  of  the  aircraft  are  symmetrically  located  about  the 
vertical  and  wing  axes  of  the  aircraft. 

3.  The  point  of  closest  approach  of  the  missile  to  the  aircraft  is  defined 
to  be  the  warhead  detonation  point. 

4.  Warhead  detonation  inside  a volume  defined  by  the  aircraft  dimen- 
sions will  have  a probability  of  hit  (Pjj)  equal  to  one. 

5.  If  the  missile  is  a hit-to-kill  missile,  a detonation  outside  this 
volume  will  have  a P^  = 0. 

6.  For  proximity  fused  missiles,  a warhead  lethality  zone  around  the 
aircraft  volume  will  be  assumed. 

7.  A warhead  (proximity  fused)  detonation  outside  this  lethality  zone 
will  have  a PH  = 0.  A detonation  between  the  two  zones  will  have  PH 
linearly  proportional  to  the  detonation  point  distance  from  the  aircraft 

volume. 

The  missile,  target,  and  flare  simulation  program  provides  miss  distance 
information  in  inertial  coordinates,  therefore,  it  is  necessary  to  perform  a 
coordinate  transformation  to  obtain  the  miss  distance  in  terms  of  aircraft 
coordinates.  Figure  8-1  shows  the  equations  used  to  perform  this  transfor- 
mation with  the  coordinates  of  the  miss  vector  being  (XMISS,  YMISS,  ZMISS) 
in  the  inertial  system  and  (XRT,  YRT,  ZRT)  in  the  aircraft  system. 
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Figure  8-1.  Coordinate 

transformation,  inertial- 
to -aircraft  coordinate 

The  aircraft  coordinate  ayatem  haa  the  X^.  and  Zj  axes  along  the 
longitudinal  and  vertical  axes  of  the  aircraft  and  Xj  axis  along  the  aircraft 
wing.  This  coordinate  system  shown  in  Figure  8-2  has  the  tailpipe  at  its 
center  and  the  aircraft  dimensions  defined  relative  to  this  point. 

For  proximity  fused  missiles,  a warhead  lethality  zone  around  the  air- 
craft volume  is  assumed.  This  zone  is  simply  determined  by  adding  to  each 
aircraft  dimension  the  warhead's  effective  kill  radius  (MR). 

If  warhead  detonation  occurs  within  the  aircraft  volume,  = 1 is 
assumed;  if  it  lies  outside  the  warhead  lethality  zone,  = 0 is  assumed.  If 
warhead  detonation  lies  between  the  two  zones,  P^  is  assumed  to  be  linearly 
proportional  to  the  distance  from  the  outer  boundary  of  the  aircraft  volume  to 
the  detonation  point.  For  the  case  of  hit-to-kill  missiles,  the  missile  effec- 
tive kill  radius  (MR)  is  set  equal  to  zero  making  the  aircraft  volume  and  war- 
head lethality  zone  coincident. 

The  equations  and  logic  required  to  implement  this  calculation  are  as 
follows: 

If  (XN  + Mr)  SXRT  and  XRT  S - (-Xfi  + MR) 

Then,  P^  = 0 

8-2 


\ 


If 


Figure  8-2.  Aircraft  coordinate  ayatom 
-(Xfi)  S XRT<  XN 


Then, 


Px'1 


Then, 


XRT  Z Xj, 


PX  = 


(XN  + Mr)  - XRT 


Mr 


Otherwise, 


PX  = 


(XB  + Mr)  + XRT 


If 


ABS  (YRT)  S (Ys  + Mr) 


Then, 

If 


PY  = 0 


ABS  (YRT)  a Yg 
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Then, 


Otherwise, 

<YS  + MR)  - ABS  (YRT) 
PY‘  Mr 

If 

ABS  (ZRT)  2 (Zs  + Mr) 

Then, 

pz  = o 

If 

ABS  (ZRT)  S Z_ 

Then, 

pzsl 

Otherwise, 

(ZS  + Mr)  - ABS  (ZRT) 

Pz  mr 

Finally, 

tr 

PH  ■ PX  ; PY  * Pz 

Figure  8-3  shows  a block  diagram  of  these  computations.  Table  8-1  gives 
the  probability  of  hit  subroutine. 


8-3.  Probability  of  hit  computation 


cr< 


Table  8-1.  Probability  of  hit  eubrouttne 
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•*•1. 

60  TO  t 
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9.  SPECTRAL  INTEGRATOR 


4 

* 

In  order  to  utilize  the  I/R  target  signature  generated  by  ASDIR  II  (or 
any  other  computer  program  which  can  generate  apparent  spectral  radiant 
intensity  1^)  It  is  necessary  to  integrate  over  the  optical  waveband  of 
the  missile  simulated  in  the  M/T/CM  program.  This  Integration  is  accom- 
plished by  use  of  an  auxiliary  routine  SPRINT,  which  can  integrate  over  any 
desired  spectral  region  of  the  ASDIR-II  output  (maximum  of  50). 

Two  options  are  available  in  SPRINT  which  can  be  selected  at  run  time 
and  extend  the  •-.-eruineas  of  this  routine.  These  include: 

(1)  An  atmospheric  transmission  table  can  be  read  in  at  run  time  and 
used  in  the  integration 

(2)  A spectral  filter  can  be  read  in  and  applied  to  the  integration. 

Using  the  optional  atmospheric  table  and/or  filter  table,  the  integrated 
JXTX  can  be  determined  as  a function  of  range  or,  if  these  tables  are  omitted, 
the  values  of  Jx'Tx**  a *unctlon  of  ran8®  fr°m  the  ASDIR-II  output  are  used. 
Figure  9-1  is  a block  diagram  of  the  SPRINT  program  and  Table  9-1  is  a 

program  listing. 


*Stone  C W..  Capt.  USAF  and  Tate,  Stanley,  Planning  ASDIR-II  (VolsI,  II, 

III)  Deputy  for  Development,  Aeronautical  Systems  Division,  Wright  Patterson 
Air  Force  Base,  Ohio,  ASD/XR-TR-75-1,  January  1975. 
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Table  9-1 


SPKINT  program  listing 


*• 

CSNH9N  SrSN(50U  SETS (50 ) < TAI5 (90 >#  HA(470)a  TBAN(470)a  *TMI0(3)a 

0010 

2. 

■ ThT(47(j)a  *XI3o6)a  WLAISOIa  WLB  (So)  a TEE  ( 50 ) a A I RX  l 50 ) a 

0020 

3. 

• *ATAK  (470 ) a TEN  TAB  (.,701  a 0SR(i70).  0S»L<1C„)a  OSSTIIOO). 

0030 

• • 

■ ST8(300 ) a M T Q ( 10a 300 1 a TAB ( 10.300 )a  anaN(A)a  LABEL ( 10*20 ) a 

0090 

s* 

• NPBEU(1O)aNa0(10)aETAB(10a300I 

r 

300 

7. 

C INITIALIZE  data 

310 

»• 

c 

320 

9. 

INTEGE9  JNITY.50URCE 

0325 

10* 

R£H|NO  9|  BE* 1*0  9(  BEHIND  7 

330 

M. 

BEHIND  10 

»• 

B1  • 3<l>159 

350 

13* 

555 

BE*0(5a30)  NSBOa ID5R  T, NPBT V, HPRTVa  XPLT VaXHJT AUaITAUH.IELTa 

1*« 

•UNITV.S9JBCE 

15. 

iE(nsho.eo.O)  stop 

U« 

OUTPUT  NSBOa IOSBT.NPBTVa ITAUH. IEIlT 

17. 

c 

340 

IS. 

C «CA0  NTAUa  LAPEL.  INPUT  BASANETS'^ 

370 

19. 

c 

310 

2P. 

IE  (UNITY. £0*1)  he  AO ( 9.5 ) NTAUi  hBIT£(3a5)  NTAui  05  T9  4 

31. 

lE(ITAUH.OT.O)  B(aD ( 5. 9 ) NT AU 

23. 

IEdTAUH.LE.OI  BEaO ( 3a 5 ) NT *U 

23. 

5 

EBBNAT (115) 

900 

2*. 

4 

CONTINUE 

25. 

OUTPUT  NTAU 

3*. 

BEAD ( 5a  1C ) (ANAhlll.  1.1.4) 

910 

27. 

10 

ESHNAT ( ?0»* ) 

*20 

n* 

wBITE ( 4. 20 ) (ANAN(I)a  I*1a4) 

*30 

29. 

20 

EOBnaT ( 1M1/32X. INON.EOUISPACEO  SPECTRAL  INTEQHaTOH i //5»Xa 'CA«E i/41 

30. 

■X.4A4///I 

950 

31. 

lE(ITAUH.ST.O)  BEAO (31  NEREQ ( 1 ) t BEA0I3)  (HA ( I ) , TRAN ( I ) , I ■ 1. 

32. 

• NEBEQIDII  00  TO  25 

33. 

25 

CONTINUE 

3*. 

NHL  .NrREQdl 

35. 

HRITEdO.lOl  ( AN  AN  ( I ) . 1.1.4) 

230 

3*. 

HB1TEI 10.30)  KPLTY 

5*0 

37. 

C 

550 

3«. 

C HEAD  INTEGRATION  LIMITS  (5  PAIBS/CABO) 

540 

39. 

C 

570 

*0* 

BEAD ( 5 a 5 0 ) (HLA(I).  hl5(I)a  IMaNSRO) 

550 

«1. 

50 

EOBNAT ( 1QE4.0) 

690 

*2* 

C 

900 

93« 

C BEAO  TABLETS  (EtLE  7#  STORE  E9R  PLOTS  5N  FILE  10 

410 

• *. 

35 

BE *0 ( 7.EN0.34  ) XT,  (LABEL  (<T,  Da  I • 1a20).NAHEO(KT) 

•5. 

XTT  • XT 

94. 

lE(XT.OT.iO)  OUTPUTlNO.  Of  TAXQETS  ExCEEOS  10»l  STOP 

620 

97. 

lE(XT.LE.O)  OUTPUT  'NO  TABOETS'I  STOP 

*30 

9|« 

OUTPUT  XTanEBEQ(XT) 

50* 

NTLT.  NEBEQ(XT) 

51. 

09  40  I . IaNTlT 

52. 

BEAD ( 7 ) hTOIXTa I ) . TA*(XTa!)aFTAB(XTa!) 

53. 

30 

EOBNAT(2*l3l 

*70 

S*. 

IE { IELT.0T.0 ) TAB(XT.l)  . ETAB(XTaI) 

55. 

*0 

CONTINUE 

5A. 

BE  *0  ( 7 ) RAOKTj 

57. 

IF(NPBTY.NE.O)  «BI TE ( 4.209)  (LABEL ( XT . I ) a 1*1.20) 

Hi 

IE(NPBTy.NE.O)  hH|TE(4a212)(hT0(XT.I  )aTAB(XTa I )a I .i.NTLT ) 

59. 

209 

E5BNaT(2C*9) 
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(Table  9-1,  continued) 

40,  tig  FOAmaTMOMTOT  L**B0A  ,11*, AH  TOT  »30X»'T0T  LAM50A»,11X* • TOT'/CSX 
|ro(PLTY. t0*0)  09  TO  94 

4*.  JF  ( (XT.NE.l ).ANO. (IOUHCE.OT ,0)1  00TB  5* 

48.  «XtTC(10,10>  (IABCI(XT,JI,  J • 1.301 

44,  HXITt(10«5S)  NTLT* ( MT0|KT« J)«TAA|XT» J),  J a 1*NTUT) 

87.  95  F0NHAT(JS/(*£lB«9Ii 

50.  5*  CONTINUE 

59.  OB  T9  35 

70«  35  CONTINUE 

71.  XT  • KTT 

7f  IF  (BOUNCE. OT.O)KT*t 

7J.  HAIT£CA,*OI  NT AU,  XT.  NBAS,  XT 

7*.  *0  FOAMaT (.15) 

75.  «AITE(A,*A8I  (al»(I J.aLBlI),  t • l’NSXQ) 

75.  5*5  FOAMAT (lOFB* a) 

77.  C X£ AO  OETECTOA  AE5X0NBE  (0FT10NAL1 

70.  C 

79.  IF(IOSXT.Cg.O)  35  TO  50 

■g.  At* 0(5, 70)  (05.1(1), OSAT(I ),  I.J,I0SNT) 

§j,  70  FOAM AT (5(F5«3»F8«A)) 

•*.  .AEOFONIF  '//IA3X.F7.  A*25x*El5.5n 

•9*  AO  CONTINUE 

•A.  C 

■7.  I main  LOAF  bn  NTau 

•5.  ........ 

19.  00  500  J.l.NTAU 

90.  IFlUNlTT.EO.il  XCAO  (5,0  ATM  10,  M0»MS*BX»0VIH 

31,  IFIITAUX.OT.O)  QO  TO  111 

9g.  IF  (UNITY. NE'l ) 00  TO  55 

53.  N4L  • NTlT 

9*.  90  SI  I«1,NTUT 

95*  TAAN(I)  a 1*0 

95.  Bt  wA (I|  a mTOII.I) 

57.  3 FOAM AT Cl A*, At 10* 3 l 

)i,  • .X|TE ( 3, AO)  n.L, atm|0, mC,hS,SX,  OvIX 

99.  NM  a Nwt/5 

100.  JF  (N4|i  »0T . (NH.5) ) nm»nm*1 

10ta  MAal 

IOC.  M|a5 

103.  90  53  XNBal.NM 

10*.  «XJT£( 3,95)  (NA("Cl,TXAN(MCI,  MCa  "A, MO) 

105.  NA  a MA*5 

105.  53  NBaMBaS 

107*  00  TO  lit 

100.  85  XEA0(3,90)  N*L,ATMI0.HC,M9,BH,0vIX 

109*  90  F0XnaT(15«3A*.5X«FB.3,5X,F8.3#5X,F8.0#8X,FA.0> 

110*  NM  a Nau/5 

111,  |F(NWt.OT. (NMa5) I NM  a NM*l 

113*  C H£ AO  HAVELEN0TM.TXAN8.  TABLE 

1 1 A*  C 

115*  "A  a 1 

115.  mB  a 5 

117.  00  100  MMtt  a 1,NM 

Hl#  XE AO  (3,95)  (HA  (MO  ,THAN (MCI,  MCaN/.»MB) 

119*  95  F0XMAT(5(FA.a-F8.A)  J 
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(Table  9-1,  continued) 


120* 
If  l* 
122* 

123. 

124. 
120* 
lilt* 
127. 
12*. 
129. 
130* 
HI  • 
132* 
133. 
13*. 
135. 
136* 
137. 
1 3ft. 
1 J9» 
1*0* 
1*1* 
1*2* 
143* 


1**. 

1*5. 

1**. 

1*7. 


!*•• 

1*». 

150. 

151. 

15?. 


iM  i 

154. 

155. 

156. 

157. 
IS*. 


159* 

1*0. 

1*1* 

162* 

163* 

15*. 


165* 

166* 

1*7. 

!*•• 

1*9* 

170. 

171. 

172. 
171. 
17*. 
175. 
17*. 
177. 
175. 
1*9* 


100 

111 

110 


"A  ■ MA.ft 
«0  • *9.5 


CONTINUE 

|F<NPNT*.NE«0)  "“ITEIAjIIO)  Ji  H*(I)«  TftANJll#  1 . 1 ,NkL  ) 

FORMAT) 1M1/I 3, 3*X# ' I NPy T MAVELENOTN  VS.  T***SM I SSI BN  •//(*(FlO«*. 


•• (U) • #C1**S> » I 
MfttTEI 10<*0I  KT.NWL 
IF(in8KT.E0*0)  o*  TO  130 
09  120  J«1.N«L 

120  9SR(|)  • TUU2(VA(1  J.OS.I  ;OS«T) 

130  CONTINUE 

c 

b COOP  ON  NUMBER  OF  TANOCTS 

c 


* • J 

tFiseuftcr «ot.oik«i 

NTUT  • NFR£U(K) 

WRITE(IO.IO)  (LABEL!*.  II.  1*1.20) 

OUTPUT  WT3|K.I).mA(1I,mT3(K.NTlT).mA (NHL ) 

OUTPUT  — 

|F(UT3(K#1).LT.mA(ii)  OUTPUTtLO*  END  OF  TUT  LESS  THAN  LAMSUA(l)'j 
•30  Tft  300 

|F ( VTQ (K.NTLT ) .UT .HA ( NmL ) ) OUTPUT  *H|OH  ENU  OF  TOT  OREaTE*  THAN  1AM 
•BOAINMU'I  STOP 
|F  (UNITY. NCI)  00  TO  135 
00  132  l.l.NTLT 
132  TtfT(l)  ■ TAN(K.I) 

K3*1|K4*NTLT 
00  TO  155 
135  CONTINUE 

00  1*0  l*l.NTLT 
STft(l)  ■ MTQIK.II 
1*0  TX  ( I ) ■ TAR(K.I) 

DO  130  I.I.NNL 

ISO  TmT ( I ) ■ TbU2(N*< I J.STO.TI) 

HI  • NTO (*. 1 ) I "2  • HTO(X.NTLT) 

CALL  SCR|PT(U1,h2.K3.K*«mA,NwL) 

155  CA'„l  TRAR ( K3.  K4.hA.ThT.  Sum  1 
AJHAX  • RAD(K)/SU" 

170  FOrJaT?///OOX#  'FRACTION  OF  TAHOET  PLUME  IN  0ANDI///50X. *HAVELEN0TM 

• RE0|0N'#7X.  »F  JTQ*/) 

00  ISO  I*1*NSR0 

Hi  • HL*< I ) . 

M2  • MLO(!) 

CALL  0C»|PT|Ml»M2,KA.XQ,MA,NMb) 

CALL  TRAP(KA,KS.MA#THT.3UMT) 

SFTS( 1 1 • BUNT/SUN 
vRITEI».175)M*<KA),w*(KS).9FTS(?} 

175  FORMAT | 50*. *3* *»T0* *«F*.J. *X.El2<5) 

ISO  CONTINUE 

30  190  1*1. NbL 
AJMAX  • l* 

190  TEMTAR(I)  • TMT < I ) *A JMAX 
00  200  Id.NSRO 
HI  • ML* (I  II  M2  • WlB(|) 

CALL  ScR|RT(Ml,U2,KA.K5.HA,NMb) 

CALLTRAP(XA.XB»HA,TEMTAR.BUMF  J 
SFSN(I)  *3UMF 

200  IF(SFSMMI.LC.O)  SFSM(I)  • 1. 


990 

1000 


1010 

1020 

1030 

1030 

1031 

1032 

1033 
103* 
10*0 
1050 
10*0 


10SO 

1090 


1100 

1110 

1120 

1130 

1131 

1132 


113* 

1135 

1139 

11*0 

1150 

11*0 

1170 

liso 

1190 

1200 


1290 

1300 

1320 

12*0 

1250 

12*0 

1270 

12S0 

1330 

13*0 

1390 


1320 

1390 

1*00 

1*10 

1*10 


(Table  9-1.  continued) 


ISO* 

1S1* 

1S2* 

710 

1S1* 

IS** 

219 

i«5* 

105* 

220 

1S7* 

ISO. 

* 

1S9. 

210 

190* 

191* 

192. 

191. 

215 

19*. 

195* 

19*. 

2*0 

197* 

19S* 

2«* 

199. 

2*5 

200* 

250 

201* 

202. 

201* 

20«* 

205. 

20*. 

207. 

201. 

209. 

210* 

211. 

75* 

212. 

213* 

21** 

755 

215* 

2*0 

215* 

217. 

2*5 

2tS* 

100 

219. 

220* 

ISO 

221* 

soo 

222* 

222. 

oa  sio  i.i.nml 

Tf»T«MU  • • THWIII 

|F (IOSRT.QT.0)  TIHTM|1»  a HHTA*  1 1 ) aDSM  1 1 ) 
m*!T£|«.719> 
raeiATci  x) 

m2|Tei*.720)  ATito.NO.HS.ie.ovia 

ra*HAT|2?l.lAA.7X.  'HO  a '.Flil.l  N'.SX.'mS  a ».M  .2.'  Hi.BX.lS*  a 
a'.M.O.'  I'.SX.'VH  a t.F*.0«'  «*'/✓) 
t>HI  TC  ( 4»?3Q  ) •AOUt.MTOlX.l  t.MTOtK.ITLT).  (L*BtL(«« I )« I al«20)»*JHAX 
*8HH*T(1«X.'KA0  a i.Mt.a.t  IN  TNC  mAVELENOTH  XEOION • .» S •«. • TO'. 
aMii.l  HIC*BNS'//1SX.iTA*5ET  TYRE  IS  '.20**//lSX. ' jSCAL  a '.E1Z.*/ 
• //> 

true.EO.o.)  veiTtis. 235)i  oa  re  2aa 

raiHAT|/t2x. •oAvrucNuTH  •taien.iox.'o  taui.iox.itau  e»f'»iox»'tau 
a (TBTAU'.YX^IAX. M -ICMSNS ' // > 
m*I TE  t (. 2*0) 

FaH1AT|/i2X. IMAVELCNOTM  *EOICN *. 10X. I J TAU ' . 10*. ' TAU  £77 i .10*. • TAU 
• |TaTAU'.9X.MHMA0tANCE'/ltX.SN|1tCeeNS)//) 

WR|TEC10.2*5)  (V*t I ). TENT Ax  1 1 ).  t a 1.Nal> 

»8MiAT<SEl*»*) 

ceiTiNur 

oa  2 AO  lal.NSNO 

Mia  ML A (III  M2aMLS<I) 

CALL  BCR ! NT ( Ml. m2. XA.xI.mA.NnL ) 

CALL  T*  AN  ( XA.KS.MA. TElTAX. SUN ) 

AATANH)  a Sul 
TEN  Cl)  a SU1/SFS"!!) 

TAleil)  a T17|I)  a S7TSCI) 

M a SMIOO. 

|F|N(.E0.0)MHlTE|«.25S)kA|KA|.MA|KS).*«TAN||).TEF|t). 
aTAISlDlOa  TO  2A0 

ramATii?x.FA,i.*.TB.i.rs.3«7x.Ei2.s.7x.rs.7.iox.rs,7) 

AlNRlJI  a AATANl  D/NS/NS 

-R!TE|*.755)MA|«Aj.aA,«B).AATA*(|).TEr|!i»TA|»|l ),Al**|j) 

FaRNAT|l2X.FS.].iaTS-l«r*.S.7X.El2>S.7X.F9.7.10X.r9>7.7X.Cl2.S| 

CONTINUE 

IFCKNJTAJ.EQ*!)  M*1TE(*.2*B)  I MAJ I ). TElTAN ( I ). lal.NWL) 
raRHATiiMl.ASX.lANLANBOA  vs.  J TAU  //(Stll.B)) 

CONTINUE 

WNtTECA.JSO) 

raRNAT(iHi) 

CONTINUE 
•a  TO  BBS 
END 


1**0 

1*70 

1AS0 

1SC0 

1510 

1520 

1S10 

1550 

15*0 

1550 

mo 

1590 

1500 

1510 

1520 

1510 

1550 

1**0 

1570 

15S0 


15*0 

1*10 

1720 

1710 

17*0 

17*0 

1770 

17*0 

1790 
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Table  9-1  (concluded) 


1' 

t< 

3. 

• t 

»• 

• • 
T# 
M 
9* 

10* 

It* 


FUNCTION  TkU?<A*K«'  } 
oinensibn  xmiiMii 
IF (XtLTtXIl ) I rLUJtFlll{«ETU»N 
09  1 I»l#50000 

IF  ( x ( I ) «3T  t X ( 1*1 ) ) TLUt»FC!>l*CTU*«N 
|FI»tOEtX(I>tANDtAtLCtXClM>»  OB  T*  * 
1 COSTIVE 

t v.(».x(i))/cxci*i)*x«in 

TUUF«F I !!**• IF  1 1*1 1-F  J I > ) 

•ETU*N 
| NO 


1* 

Ft 

3t 

• t 
»t 
»t 
?• 

• t 

9* 

10t 

u» 

IF* 

S3* 

!•# 

I9t 


9U*MuTInE  •Ci*IFT|ni##J#LI#LJ#navE#n*L) 
OINJNSION  -AVtlll 
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10.  SAMPI.F  RUNS 


A act  of  runs  showing  the  interaction  of  the  entire  methodology 
simulation  (ASDIR  II,  M/T/CM,  end  SPKINT > was  made,  and  sample  outputs 
are  given  in  this  section.  The  engli.  • osed  to  calculate  the  IR  signature  was 
the  10,  000  foot  default  engine  of  ASOIR  II.  Engine  hot  part  contributions 
were  assumed  using  the  equivalent  blackbody  temperature  and  area  of  824°K 

. 2 

and  730  cm  respectively  (0  degrees  aspect  angle).  The  blackbody  area  was 
varied  as  a function  of  the  cosine  of  the  aspect  angle  which  was  varied  at 
15  degree  increments  from  9 to  90  degrees.  Apparent  values  (1.8  to 
5.5  microns)  were  calculated  for  ranges  of  0.0,  0.305,  A. 524,  6.096,  and 
15. 240  kilometers.  These  values  were  then  integrated,  by  SPKINT,  over 
five  spectral  intervals  to  generate  the  apparent  effective  (Jt)^  values  used 
in  the  M/T/CM  simulation  program.  The  integrated  values  of  (JT)^  were 
then  entered  into  the  M/T/CM  program  and  a typical  set  of  missile  simula- 
tion runs  using  a spin-scan  type  missile  were  made  for  aspect  eagles  in 
15  degree  increments  from  8 to  90  degrees,  and  launch  range  of  5,  000  feet. 
The  results  of  these  runs  were  = 0 for  all  but  aspect  angleo  of  75  and 
90  degrees.  At  these  two  angles,  the  missile  was  unable  to  maneuver  to 
catch  the  target  (i.e.,  it  was  launched  outside  of  the  aerodynamic  launch 
boundary)  and  P#  * 1 in  these  cases. 

Sample  outputs  from  these  runs  are  shown  in  Tables  10-1  and  10-2  and 
in  Figures  10-1  through  10-6.*  Table  10-1  shows  the  ASDIR  II  output  of 
\TX  versu*  ^ Table  1-2  the  integrated  SPKINT  values.  Both  of  these 
cases  are  for  0 degrees  aspect  and  0 Km  range.  Figure  10-1  gives  a plot  of 
the  spectral  J^-r^  for  the  0 Km  range  case,  and  Figure  10-2  is  a polar  plot 
of  ^ ~ 1.8  to  2.  6p,  for  three  ranges.  Plots  of  the  simulated 

missile  flight  are  shown  in  the  last  four  figures.  Missile  target  trajectories 
in  the  X-Z  and  X-Y  planes  are  shown  in  Figures  10-3  and  10-4  respectively. 
Apparent  effective  Intensity  and  effective  irradiance  at  the  missile  seeker  as 
a function  of  time  are  given  in  Figures  10-5  and  10-6. 


*Figures  10-1  through  10-6  were  generated  by  separate  CALCOMP  plotter 
routines  which  are  not  part  of  EPICS. 
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Table  10-1.  ASD1R  II  output 
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Table  10-2.  SPKINT  integrated  output 


•me 


*0 


lilOfta*  ;>t)t 
2*672»T8*  J.bSO 
2*S»1»T8*  *.777 
3«S32»T8*  *.888 
3.I58-T8*  b*  1*8 


3*05  M 


Mb  ■ 


. 18858E 
•31717C 
•818S5E 


00 

00 

00 


•20835E  00 
•38818E  00 


•3. OS  " 


•8  • 


Q.  1 


V*  • 


o. 


8*0  • .28131  03  Ifc  TM£  W*vElIn0Tm  8£oJ8N  1*80*0  T8  5**25*  «IC"8N» 

TA80CT  T <lPt  IS  *83|8  TEST  8(J*  10«*T  0£7*UUT 
JSC*U  • *10001  01 


«*veuen3tm  «tans 
MIC88N8 


J TAJ 


TAU  1.7* 


T*U  ( TOTAL 


1.808*T8*  2*838 
2*87?«T8«  3*5*0 
2*88fTb.  *.777 
3*S3?.T8-  **888 
3*858*T8»  5*1*8 


02 

02 


•*87*7E 
•87510E 
•17o8*E  03 
•S7*8*E  02 
•10737E  03 


1.0000000 

1*0000000 

1*0000000 

1*0000000 

1*0000000 


*1888770 

•3171721 

•8135S03 

•208352* 

•3881558 


10-3 


- '-.*•  ' * » *"  -Mfe.  ».  • 


Figure  10-1.  Spectral  (R  = 0,  aspect 
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Figure  10-2.  PoUr  plots  of  apparent  effective  radiant  intensity 
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Figure  10-5.  Apparent  radiant  intensity  at  miss'le  seeker 
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